TABELA 2 


Massas especificas de soluções aquosas de álcool 


x (9%) 


PexpíO em?) 
PoatÃO em *)* 


(4p/p) 100 ** 


[09878 | 09847 


0,98776 | 0,98475 | 0,9819 


E 0,004 ess 00, 


* Nos valores de Peste retivemos um algarismo a mais por 


se tratar de valor a utilizar em cálculos subsequentes. 


ii dp = Pexp “Prale 


Consideremos, a titulo de exemplo, uma série de 
determinações da massa especifica de misturas de água 
e álcool, a 20' C, cujos resultados estão contidos na 
Tabela 2. A composição das misturas é representada pela 
fracção em massa, x, do álcool. 

O valor correspondente a x = O merece particular 
confiança, pois a massa especifica da água pura a 20 
é conhecida com uma exactidão da ordem de 1/10º. Por 
outro lado os valores de p para x >0 são relativos a 
esse padrão (água pura) e as condições experimentais 
permitiram apenas assegurar uma reprodutibilidade da 
ordem de 1/10. 

O gráfico da fig. 1 apresenta uma ligeira curvatura, 
pelo que a equação 


p=Pf+ ax (8) 


não é a mais apropriada para traduzir os resultados 
experimentais. 


FIG. 1 


Gráfico de p em função de x. A tracejado, a «recta que melhor 
se ajusta» 
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Mas (fig. 2) o gráfico de (p; — p,)l(x, — x5) em fun- 
ção de (x;+xjé sensivelmente linear; podemos, pois, 
tentar a equação 


p 
l (e; = Poll; fa X,) e 


Gta, x+ ax (9) 


Fazendo F 


Il 


G, = X FT Xy 


o problema consiste em determinar os parâmetros da 
equação F = a, + a;G. Tabelando os valores de F, e 


e 
Eai O DR od > + >> + — e 


2 4 5 8 O = 


Gráfico de (ppl) em função de x +x, 


de G, e utilizando as fórmulas (1), (2) e (7) calculamos 
os parâmetros a, = — 1,905 X 10º gem? 


a=2,75x10ºg em? e p= 0998209 em”. 


A consulta da Tabela 2 permite verificar que a 
equação (9) reproduz muito satisfatoriamente os resul- 
tados da experiência, pois os desvios relativos não 
excedem 0,004% (fig. 3) e o desvio relativo médio é 
0,001; 9%. 

A equação (8) com p,=0,99783g cm“ e a= —1,63X 
x 10!g em, também obtidos por análise de regressão 
linear, só deve ser adoptada quando não for necessário 
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grande rigor. Os resultados experimentais são reprocu- 
zidos com desvios relativos que podem atingir 0,04 % 
(fig. 3) e o desvio relativo médio é 0,025 %, ou seja, 
vinte vezes maior que o correspondente à equação (9). 


GO 


x 
g02 
G 
E 
“= e 
bo) a Rai o e ' ã é + a 
EN a [=] Mm 
da af e FE 
Mao x x 
-BoZ 
us x 


FIG. 3 


Gráfico dos desvios relativos. X cg 8 “cega 

Os cálculos referentes a este problema foram reali- 
zados em poucos minutos com uma calculadora portátil 
equipada para realizar análises de regressão linear. Al- 
gumas calculadoras portáteis estão equipadas para for- 
necer também o coeficiente de correlação entre x e y, 


mas os desvios padrão dos parâmetros calculados só 
são facilmente acessíveis se se dispõe de meios algo 
sofisticados de cálculo: calculadoras de secretária de 
certo porte ou computadores. Contudo, estes requisitos 
são dispensáveis em muitos casos, pois o desvio rela- 
tivo médio e um gráfico de desvios relativos dão uma 
imagem bastante fiel da dispersão dos valores experi- 
mentais relativamente aos calculados a partir da equação 
adoptada. 
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TRABALHOS 
DE REPARAÇÃO 
E CONSOLIDAÇÃO 
DA BARRAGEM 
DO ROXO 


E 
e 


Tratamento do talude de contenção da albufeira, 
no encontro da margem direita, do lado montante, 
de forma a permitir a escavação com cerca de 
23 m. de profundidade para o prolongamento da 
barragem gravidade. 


Natureza dos trabalhos realizados: 


e Cortina plástica de estanqueidade, moldada no 
solo. 

e Sondagens — amostragem integral. 

e Furação para injecções — ensaios tipo Lugeon. 

º Injecções para impermeabilização do maciço 
rochoso no prolongamento da cortina plástica. 

e Instalação de piezómetros e inclinómetros. 

e Drenagem de maciço rochoso. 


TEIXEIRA DUARTE, LDA. 


SEDE: 


Av. da República, 42 / Telef. 73 30 86 / Telex 18548 TEDEX P / 1094 LISBOA Codex 
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O consórcio Brown Boveri 
na distribuição de energia 


O disjuntor é um O nosso programa de Planificamos e cons- 
elemento base na fornecimento compre- | truimos instalações 
distribuição de energia  ende: exteriores e interiores 
eléctrica — Disjuntores pneu- para todas as tensões 
, máticos ate 765 kV 


— Disjuntores de peque- 
no volume de óleo 
— Disjuntores SFs 


Os disjuntores BBC Condicionamentos Com prazer resolve- 
respondem a todas técnicos, de serviço e mos o vosso problema 
as exigências de econômicos estão na 

serviço base da escolha de um 


dos três sistemas de 
extinção do arco 
eléctrico 


1531564 


> o Ss » 
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o ASA EA a a ea > e Sie e 13 $ 
2 LOST: A: Perde Ds poses ra ATA ” Za na o 
15.4 124% 44 Air, 2a trato £ f 
ES ISTO TA A ASAE USA A POD O SAIAS SD ADÃO PSA E 


Disjuntores pneumá- 
ticos 


Disjuntores de Disjuntores SF, 
pequeno volume de 
óléo 
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Ocorrência média anual no território Português 
das classes de estabilidade atmosférica Pasquill-Giflord 


RESUMO 


Com o objectivo de fundamentar a regulamentação 
aplicável às actividades que contribuem para a poluição 
atmosférica no Território Nacional, efectuou-se a caracte- 
rização da ocorrência média anual das classes de esta- 
bilidade atmosférica Pasquill-Gifford no Território Conti- 
nental, 

Esta caracterização permitirá estabelecer as expres- 
sões mais adequadas para o cálculo da dispersão de 
efluentes atmosféricos no Pais. 


OBJECTIVO 


— Fundamentar a regulamentação aplicável às acti- 
vidades que contribuem para a poluição atmosférica no 
território nacional. 

— Pela sua natureza genérica as medidas a promulgar 
devem abranger actividades poluidoras indiscriminadas de 
pequena dimensão e baixa perigosidade quanto aos eflen- 
tes que emitem. A regulamentação deve ser simples de 
aplicar e terá como principal objectivo evitar que sejam 
excedidas concentrações médias anuais admissiveis. São 
especificamente excluídas unidades poluidoras de grande 
porte ou de elevada perigosidade, 

— Na concretização do objectivo proposto é finali- 
dade deste trabalho a caracterização da ocorrência média 
anual das classes de estabilidade atmosférica de acordo 
com a tipologia internacionalmente mais utilizada (Pas- 
quill-Gifford), de modo a permitir estabelecer as expres- 
sões mais adequadas para o cálculo da dispersão de 
efluentes atmosféricos no território português, 


SUMÁRIO 


Determina-se a frequência média de ocorrência das 
classes de estabilidade de Pasquill-Gifford para o terri- 
tório português com base em observações de superficie 
recolhidas pelo Instituto Nacional de Meteorologia e Geo- 
fisica nos anos de 1970, 1971, 1972 e 1973. 


J. J. DELGADO DOMINGOS 
M. FATIMA PINTO 
M. TERESA PONTES 


ABSTRACT 


Annual mean occurence of Pasquill-Gifford stability 
classes in Portugal was made in order to fundament the 
policy to rule atmospheric pollutant activities, 

Based on this caracterization more correct expres- 
sions to calculate the dispersion of atmospheric pullutants 
in Portugal can be set, 


Estima-se que os desvios médios anuais não ultra- 
passem em mais de 5% os obtidos ao considerar um 
número de anos mais elevado, 

Tal precisão está de acordo com a finalidade do 
trabalho e a precisão máxima esperável na estimativa da 
dispersão de poluentes usando a metodologia internacio- 
nalmente mais difundida. 

A tipologia subjectiva de Pasquill-Gifford foi substan- 
cialmente melhorada através da caracterização objectiva 
de Turner cuja aplicação se discute, 

O trabalho considerou cerca de 49 milhões de dados 
recolhidos em 160 estações meteorológicas do |. N.M. G. 
Para cada estação apresentam-se, ano a ano, os valores 
médios trimestrais e anuais. 

Com base nos valores obtidos, estabelecem-se gran- 
des zonas regionais de ocorrência média anual das dife- 
rentes classes, discutindo-se a validade e aplicação das 
conclusões obtidas. 


1, INTRODUÇÃO 


11 — O estabelecimento de bases adequadas para 
a previsão e controlo da poluição atmosférica exige o 
conhecimento das propriedades de transporte e dispersão 
da atmosfera relativamente aos poluentes que nela são 
lançados. 


Trabalho realizado para o Serviço de Estudos do Ambiente da Secretaria de Estado do Ambiente — 1978. 


Manuscrito recebido para publicação em 31/3/B0. 
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As características de transporte e dispersão da atmos- 
tera são também determinantes em fenómenos vitais da 
biosfera, pelo que um planeamento biofísico fundamen- 
tado deve, necessariamente, tê-las em conta. Efectiva- 
mente, a atmosfera constitui um elo fundamental nos 
ciclos biogeoquimicos naturais (ciclo da água, do carbono, 
do azoto, do enxofre, etc.), cujo elemento motor é o Sol. 
Na conversão da energia solar em correntes de trans- 
porte, a atmosfera desempenha papel análogo ao do fluido 
termodinâmico num motor térmico. 


Para que se ultrapasse uma visão meramente quali- 
tativa relativamente aos ciclos fundamentais, é essencial 
saber caracterizá-los não apenas em relação às quanti- 
dades totais envolvidas mas sobretudo quanto aos fluxos 
(massas por unidade de tempo e superfcie) que se 
verificam a cada instante, A poluição atmosférica é, 
actualmente e na maioria dos casos, a consequência de 
o fluxo de um composto determinado exceder a capaci- 
dade que a atmosfera tem de o dispersar a velocidade 
suficiente por todo o planeta de modo a não serem 
excedidas as concentrações que ultrapassam a capaci- 
dade regeneradora local da biosfera. 


t.2— A dispersão de espécies quimicas gasosas, de 
humidade, de calor, de particulas em suspensão, etc., é 
basicamente determinada pela turbulência da baixa atmos- 
fera. À turbulência consiste na flutuação aleatória da 
velocidade de um fluido em torno do valor médio obtido 
num dado intervalo de tempo, Esta velocidade média tem 
as caracteristicas de um transporte à qual se sobrepõe 
a dispersão que promove a sua diluição. 

A turbulência é um fenómeno dissipador de energia 
mecânica pelo que a sua permanência exige uma fonte 
positiva de energia. Tal fonte é, em última análise, o Sol. 
Pode, todavia, considerar-se um efeito directo, resultante 
das variações de temperatura do solo devido às per- 
mutas de energia radiante com o Sol e o ambiente 
circundante, e um efeito indirecto originado na própria 
existência do vento e da interacção deste com a super- 
ficie terrestre. 


O valor relativo da contribuição para a turbulência 
dos efeitos mecânicos do vento e dos que resultam de 
variação de temperatura da superfície e suas diferenças 
com o ambiente exterior, determina o que é habitual 
designar por estabilidade da atmosfera. À dominância dos 
fenômenos térmicos na geração da turbulência dá origem 
a atmosferas muito instáveis. À dominância de efeitos 
mecânicos, atmosferas neutras, A supressão ou amorte- 
cimento da turbulência devido a variações inversas de 
temperatura origina as atmosferas estáveis ou muito está- 
veis. À atmosfera muito instável corresponde o máximo 
de dispersão local, À muito estável, o minimo. 


2. ESCALAS ESPACIAIS E TEMPORAIS 
NO TRANSPORTE E DISPERSÃO DE POLUENTES 


21 — À turbulência verifica-se no escoamento de 
qualquer fluido a partir de certos valores da sua veloci- 
dade e corresponde a um fenómeno de instabilidade de 
caractersticas não lineares. É um fenómeno extremamente 
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complexo, objecto de investigação teórica e experimental 
desde há dezenas de anos. À sua complexidade resulta, 
em grande parte, das suas características não lineares, 
as quais provocam interacções entre as flutuações da 
velocidade do fluido (e das espécies quimicas que trans- 
porta) envolvendo escalas espaciais que vão da fracção 
de micron a milhares de quilômetros, e de fracções de 
microsegundo às muitas centenas de horas quando se 
trata da atmosfera. 


No âmbito deste espectro, quando o dominio de 
de aplicação é a atmosfera, a gama até 1—2Km é 
usualmente considerada como a da micro-meteorologia; 
até — 100 Km, a da meso-escala; para além de — 100 Km 
a da macro-escala ou escala sinóptica. 


22 -—No estudo da poluição atmosférica, é neces- 
sário ter em conta todas as escalas quando se trata de 
abordar a evolução global de um poluente, Tal análise 
global raramente é abordada, muito embora seja crescente 
o seu interesse, nomeadamente quanto ao estabelecimento 
da responsabilidade que uns países têm na poluição que 
noutros se verifica. O estudo dos fenómenos de meso- 
-escala é importante para o próprio pais, não apenas 
sob o aspecto de transporte inter-regional de poluentes 
mas também porque tais estudos base são de aplicabili- 
dade imediata num planeamento biofísico com base 
quantitativa. 


As escalas referidas, embora interdependentes, têm 
características próprias que permitem separá-las numa 
primeira abordagem desde que as limitações que de tal 
decorrem sejam claramente entendidas. 


2.3 — Às escalas espaciais estão associadas escalas 
de tempo que se podem grosseiramente ligar ao tempo 
que medeia entre a ocorrência de um fenómeno num 
local e a sua detecção noutro, Tal escala temporal é da 
ordem das dezenas de horas à escala sinóptica facto que 
se encontra na base da previsão meteorológica corrente 
nessa escala tirando partido do facto de, através dos 
modernos sistemas de telecomunicação ser possivel, 
de modo praticamente instantâneo, transmitir a informação 
relevante obtida nos postos de observação que cobrem 
a escala espacial correspondente. À escala sinóptica 
estão associados factores físicos dessa escala e são 
esses que a previsão meteorológica corrente prevê com 
tanto mais rigor quanto menor for a escala de tempo 
abrangida pela previsão. Tal deve-se ao facto da inte- 
racção da evolução sinóptica com fenômenos de meso- 
-escala se tornar cada vez mais importante até a modi- 
ficar substancialmente. Por sua vez, os fenómenos de 
meso-escala são progressivamente afectados pela sua 
interacção com a micro-escala, Esta interacção, de duplo 
sentido, entre as várias escalas é, pela sua complexi- 
dade, o factor determinante na incapacidade actual de 
previsão a longo prazo (semanas ou meses) da evolução 
meteorológica global, 


24 -—No estudo da dispersão de poluentes ao nivel 
da micro-escala e da meso-escala é fundamental o conhe- 
cimento da turbulência na baixa atmosfera na qual são 
dominantes os efeitos de superficie. Tais efeitos esten- 
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dem-se verticalmente e de modo predominante até à 
ordem do quilômetro, Embora a estrutura desta zona 
seja dominada por efeitos de superficie, a sua evolução 
global é determinada pela evolução sinóptica, Todavia, 
esta é significativamente mais lenta pelo que tal evolução 
é habitualmente considerada com um dado externo no 
estudo de efeitos locais. Acentue-se que a evolução 
sinóptica recebe também a influência de efeitos locais, 
muito embora de uma forma integrada que atenua carac- 
terísticas pontuais, Existe assim, uma jerarquia de causa 
efeito: as pequenas escalas influenciam o comportamento 
a escala global mas as suas peculiaridades locais 
diluem-se no conjunto, Às grandes escalas determinam 
a evolução (lenta) do que se passa nas pequenas escalas, 
Nestas, são dominantes os efeitos locais, 


2.5 — Tendo em conta as escalas referidas, é com- 
preensivel o facto de a meteorologia clássica se ter 
sobretudo preocupado com a escala sinóptica, só recen- 
temente começando a aprofundar o conhecimento dos 
fenômenos na camada limite atmosférica, Por outro lado, 
o estudo da turbulência e dos escoamentos em camada 
limite é de há cerca de cem anos objecto de extensa 
investigação teórica e experimental noutros dominios da 
ciência aplicada e da engenharia, nomeadamente na 
aerodinâmica, na mecânica dos fluidos reais, na transfe- 
rência de calor e massa, etc, os quais se entroncam 
no aspecto mais fundamental da Termodinâmica (Aplicada) 
dos processos irreversíveis, Às aquisições nesse dominio 
são imediatamente aplicáveis à camada limite atmosférica, 
numa interacção interdisciplinar que permitiu já avanços 
consideráveis. Nesta perspectiva, o conhecimento actual 
mostra que a caracterização das propriedades locais de 
dispersão da atmosfera não exige a proliferação de 
postos de observação meteorológica destinados a acumu- 
lar biliões de dados destinados a eventual tratamento 
estatístico. Uma análise aprofundada revela, pelo con- 
trário, que muitos postos são desnecessários e inúmeras 
medidas irrelevantes, muito embora algumas observações 
cruciais, de natureza pontual, continuem a não ser feitas 
Sob este aspecto, Portugal não é excepção, 


2.6 — Conhecida a evolução à escala sinóptica, é 
viável prever a evolução na meso-escala e na micro- 
-escala usando o conhecimento fisico relevante e uma 
adequada formulação e tratamento matemático. Nesta 
base, muitos postos de observação tornam-se desneces- 
sários, Para efeitos de planeamento, não é a previsão 
da evolução num dado momento o que se procura esta- 
belecer nas sim a caracterização estatistica de ocorrên- 
cias posteriores baseadas em conhecimento passado, Para 
esta finalidade, as evoluções sinópticas passadas são, 
em princípio, conhecidas. Não o são, todavia, situações 
locais, pois que ao nível da micro-escala é inviável dispor 
de informação estatística suficiente devido ao número de 
postos de observação que seria necessário ter mantido 
em funcionamento, Por outro lado, a instalação de um 
posto meteorológico num local onde se prevê a insta- 
lação de importante complexo poluidor, raramente será, 
por si só, capaz de efectuar, em tempo útil, a obtenção 
de informação estatística efectivamente relevante pois tal 
exige vários anos de observação. Se a informação pre- 
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tendida é apenas de nivel genérico, o presente trabalho 
mostra que o número de anos necessário é reduzido se 
apenas de valores médios anuais se trata. Nesse caso, 
porém, bastará em muitos casos o tratamento já efectuado 
no âmbito deste trabalho, 


3. AS CLASSES DE ESTABILIDADE DE PASQUILL 


3.1 — À caracterização das propriedades dispersivas 
da atmosfera internacionalmente mais vulgarizada é feita 
em termos das chamadas classes de estabilidade de 
Pasquill (em homenagem ao seu autor, que as propos 
em 1961 [1]) as quais se tornaram conhecidas por Pas- 
quill-Gifford após a reformulação por este apresentada [2]. 
Tal metodologia é utilizada em praticamente todos os 
estudos referentes à dispersão de poluentes radioactivos 
(de Centrais Nucleares, em funcionamento normal ou de 
acidente, armas nucleares, etc.), de efluentes gasosos ou 
particulas sólidas que fundamentam regulamentos oficiais 
de segurança nuclear ou de prevenção e combate à 
poluição atmosférica. 

O método foi extensamente testado com resultados 
experimentais e constitui uma base de previsão com rigor 
aceitável no âmbito da micro-escala, embora a sua exten- 
são à meso-escala seja frequente. 

Apesar das suas intrinsecas limitações e de frequen- 
tes e abusivas extrapolações na consideração de casos 
individuais, a finalidade especifica deste trabalho justifica 
a sua utilização, inclusivamente pela facilidade com que 
se estabelecerão pontes fáceis com legislação existente 
noutros países. 


32-— A estabilidade atmosférica na tipologia de 
Pasquill-Gifford abrange seis classes, que vão da muito 
instável a muito estável, Como referido na introdução, à 
classe instável corresponde a maior dispersão local e à 
estável a menor, Em termos aproximados, há uma rela- 
ção sensivelmente linear na intensidade da turbulência 
entre as classes extremas, a qual varia cerca de dez 
vezes entre os extremos. 

Na tipologia de Pasquill-Gifford apenas intervêm 
observações meteorológicas correntes. No quadro seguinte 
apresenta-se tal caracterização: 


Classes de estabilidade de Pasquill-Gifford 


A — atmosfera muito instável 

B — atmosfera moderadamente instável 
C — atmosfera ligeiramente instável 

D — atmosfera neutra 

E — atmosfera moderadamente estável 
F — atmosfera muito estável 


Velocid vento | INSOLAÇÃO DIURNA CONDIÇÕES NOCTURNAS 


à Sup. Forte Moderada Fraca Nebulosidade Nebulusidade 
(m/s) > 4/8 <3/8 
<2 A A-B B 

2 A-B B C E F 

4 B B-C GC D E 

5 Q C-D D D D 
6 C D D D D 


29 


A vantagem da tipologia de Pasquill-Gifford é a 
facilidade da sua aplicação e a simplicidade dos con- 
ceitos físicos subjacentes. Tal facto, associado às corre- 
lações empíricas que foram estabelecidas tomando-a como 
base, estão na origem da sua rápida generalização. 
Todavia, a formulação de Pasquill-Gifford é essencial- 
mente qualitativa ao considerar os factores determinantes 
das caracteristicas dispersiva da atmosfera, Nomeada- 
mente, a caracterização da permuta de energia entre o 
solo e a atmosfera é reduzida a termos subjectivos como 
insolação diurna forte, moderada, fraca e a trocas de 
energia radiante durante a noite fazendo intervir apenas 
nebulosidade de modo simplificado. Na caracterização 
dos factores de origem mecânica intervém apenas a 
velocidade do vento, ignorando-se o efeito de rugosidade 
de superficie, factores locais de relevo, etc, 

Às intrínsecas limitações referidas, acresce que a 
subjectividade inerente aos designativos de forte, mode- 
rada, fraca não permitem uma caracterização uniforme 
e objectiva mesmo no âmbito da tipologia usada, 


3.3 — De modo a resolver as ambiguidades na carac- 
terização de Pasquill-Gifford, uma definição objectiva foi 
introduzida por Turner [3], na qual a quantificação da 
permuta radiante de energia é significativamente melho- 
rada tendo em conta a energia solar incidente e a altura 
de nuvens, Muito embora permaneçam outras limitações, 
inerentes à simplicidade e generalidade do método, a 
caracterização objectiva de Turner, além da eliminação 
de aspectos subjectivos, presta-se melhor ao tratamento 
automático de grande volume de informação necessária 
à obtenção de significância estatística. 

Na caracterização de Tumer, são introduzidas sete 
classes de estabilidade, numeradas de 1 a 7. Uma análise 
aprofundada realizada por Golder [5] entre as classes 
de estabilidade de Pasquill e as de Tumer, mostrou a 
seguinte equivalência como sendo a mais adequada: 


QUADRO | 


Correspondência entre as classes de estabilidade 
de Pasquill-Gifford e de Turner 


| Pasquill-Giford Turner 


Classes de Estabilidade 


Muito instável 
Moderadamente instável 
Ligeiramente instável .. ... 
BIB us mamas dE ti ad 
Moderadamente estável... 
Muito estável 


amoo » 
Jorn 


Verifica-se assim que Turner introduz uma classe inter- 
média — 5 — que se situa entre a ligeiramente estável e 
a neutra, de Pasquill, 


3.4 — Admitida a correspondência anterior entre a 
caracterização de Pasquill e de Tumer, levanta-se obvia- 
mente a questão de saber que vantagens há em introduzir 


a tipificação de Turner, ou qual o interesse em reverter 
de uma caracterização de Turner para a de Pasquill, 

A filosofia adoptada neste trabalho parte da consta- 
tação de que tanto a tipificação de Pasquill como a de 
Tumer são insuficientes para a caracterização completa 
das propriedades dispersivas da atmosfera devido ao 
facto de nenhuma ter em conta a totalidade dos parã- 
metros relevantes. Deste modo, a uma mesma classe, 
seja qual for o tipo de caracterização adoptada, nunca 
corresponderão exactamente, em locais diferentes (nem 
sequer no mesmo local) dispersões iguais, Entre os 
valores previstos e observados existirão sempre desvios. 
Estes desvios são inerentes à simplicidade do método. 
Para a finalidade deste trabalho esses são, todavia, 
aceitáveis. Todavia, se os desvios são aceitáveis para 
uma definição das classes de Pasquill-Gifford efectuada 
por um meteorologista experiente que subjectivamente 
corrige a ausência de separações marcadas entre classes, 
já o mesmo não sucede se o volume de dados a tratar 
for muito elevado ou se a classificação é feita para 
locais muito diferentes onde não houve medidas de 
campo. Sendo objectiva, a classificação de Turner remove 
aquelas ambiguidades. Reverter da classificação de Turner 
à de Pasquill-Gifford através da correspondência esta- 
belecida por Golder tem o interesse de tornar mais fácil 
a comparação com uma tipologia largamente utilizada. 
Neste trabalho, adoptou-se a caracterização objectiva de 
Tumer, À correspondência de Golder será utilizada apenas 
para referência e ponte para a vasta literatura em que 
as correlações se exprimem em termos das classes de 
Pasquill-Gifford. Efectuar ou não a correspondência em 
casos particulares de aplicação depende da informação 
disponível localmente e do tipo de aplicação. Refere-se 
todavia que a caracterização de Turner é utilizada geral. 
mente em estudos de poluição urbana e mais recente 
mente em estudos de segurança de Centrais Nucleares, 


4. CARACTERÍSTICAS DA INFORMAÇÃO DE BASE 
E SEGURANÇA DOS VALORES MÉDIOS OBTIDOS 
PARA O TERRITÓRIO NACIONAL 


41 — A análise efectuada, de natureza estatistica, 
toma como base informação de superficie recolhida 
pelo Instituto Nacional de Meteorologia e Geofisica em 
160 estações meteorológicas, sendo 126 no Continente 
e 34 nas Ilhas Adjacentes. 

A informação tratada corresponde aos anos de 1970, 
1971, 1972 e 1973, por ser aquela para a qual havia já 
registo em suporte magnético, 

A localização dos postos de observação encontra-se 
assinalada no Mapa |. Para efeitos de tratamento e refe- 
rência posterior, cada estação é identificada pelo seu 
número mecanográfico. À correspondência entre o número 
a localidade e as suas coordenadas geográficas, constam 
do Anexo | No Anexo |* apresentam-se, estação a esta- 
ção, as coordenadas geográficas, o número de obser- 
vações, os valores mensais, trimestrais e anuais da 
frequência relativa com que occorreram as várias classes 
de estabilidade, as quais foram agrupadas de 1 a 7 de 
acordo com a caracterização objectiva de Turner, 


* O Anexo | não é apresentado dado o enorme volume de informação que contém, 
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42 — À informação referente a quatro anos corres- 
ponde a um volume total de aproximadamente 49 milhões 
de dados. Destes, apenas um subconjunto de 7,5 milhões 
é directamente relevante para o estabelecimento das 
classes de estabilidade. Todavia, antes de ser efectuado 
qualquer tratamento foi necessário efectuar sobre toda a 
informação testes prévios de consistência e validação 
de modo a detectar anomalias, A detecção de anomalias 
apola-se em critérios indirectos, os quais se dirigem 
sobretudo a erros de registo, 


4.3 — Ao efectuar-se um tratamento com a dimensão 
que facilmente se descortina, mesmo usando o cálculo 
automático, a questão central que naturalmente surge é 
a do significado final dos valores obtidos. Tal significado 
encontra-se estreitamente ligado à qualidade de dados, 
à frequência da sua recolha e ao periodo total a que 
se referem, tendo em conta o objectivo final que se 
pretende. 


Como as anomalias detectadas em alguns casos 
foram corrigidas com a colaboração do |. N. M. G., quando 
se deviam a erros de registo e o seu número não tem, 
estatisticamente, significado, a qualidade dos dados não 
é posta em causa, 


44 -— Relativamente à frequência das observações 
diárias, elas variam entre 24/dia (horárias) e 2, com casos 
intermédios de 8 e 3. Por outro lado, sendo fixas as 
horas de observação, certos valores correspondem num 
dado mês e dia a periodo diurno e noutro ao período 
nocturno, Com os dados existentes, resulta que para 
muitas estações, o número de observações anuais corres- 
pondente ao perodo diurno é significativamente superior 
ao periodo nocturno. Ao estabelecer-se a ocorrência mé- 
dia anual das várias classes de estabilidade, é necessário 
ter em conta a observação anterior. O esquema de 
correcção adoptado é descrito no capitulo seguinte, 
devendo no entanto acentuar-se, desde já, que a carac- 
terização das propriedades dispersivas das classes de 
estabilidade que ocorrem durante a noite é a que maior 
incerteza oferece devido à grande variabilidade de situa- 
ções que podem ocorrer, Todavia, tal questão não pode 
ser razoavelmente solucionada com o simples aumento 
de frequência das observações desde que se mantenha 
a tipologia de Pasquill-Gifford. Aliás, o mesmo acontece 
com os periodos de amanhecer e de anoitecer, em que 
o esquema simples de PasquilkGifford é manifestamente 
insuficiente, apesar de extensamente utilizado internacio- 
nalmente. As limitações anteriores mostram que a mera 
existência de um grande número de observações, do tipo 
usado para a caracterização de Pasquill-Gifford, não 
aumenta por si só a precisão como há a tendência 
errônea de considerar. 


Excluindo as classes nocturnas (e de um modo geral 
as muito estáveis), a precisão geral esperável numa tipo- 
logia do tipo adoptado é compativel com a menor fre- 
quência de observações nas estações não sinópticas. 
Tal deve-se a que a transição de uma classe de esta- 
bilidade a outra se processa a uma escala de tempo que 
as observações sinópticas permitem caracterizar. Nestas 
circunstâncias, ao definir por interpolação grandes zonas 
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de classes de estabilidade, a ponderação diferente com 
que figuram as estações sinópticas introduz um efeito 
correctivo que se admite adequado para a finalidade 
deste trabalho, A consulta dos valores individuais para 
cada estação (que figuram em anexo) pode, em casos 
particulares, permitir o esclarecimento de dúvidas que 
eventualmente surjam, Assinala-se, todavia, que a carac- 
terização pontual de condições locais apenas baseada 
em valores obtidos nesse ponto não aumenta necessa- 
riamente o rigor da caracterização dispersiva da «atmos- 
fera local» e pode levar a erros grosseiros sobretudo 
quando por essa via se pretendem ultrapassar as limita- 
ções intrinsecas à tipologia de Pasquill-Gifford, atribuin- 
do-lhes um rigor que de facto não possui, Nessas 
circunstâncias, a caracterização terá de adoptar esquema 
diferente, 


4.5 — Para além do rigor dos dados e da frequência 
das observações importa estabelecer a significância esta- 
tistica dos valores médios anuais, que constituem a 
finalidade deste trabalho. 


Segundo Pasquill, o cálculo da dispersão de efluentes 
em condições bem definidas de estabilidade é conside- 
rado como muito bom se em médias de longa duração 
a diferença entre valores estimados e medidos não excede 
os 10— 15%. Para médias obtidas durante intervalos 
curtos (entre poucos minutos e a hora) tal diferença pode 
situar-se na relação de 1 para 2. Estes valores acentuam, 
naturalmente, a importância do periodo médio considerado 
e estabelecem, indirectamente, um limite para a margem 
de incerteza admissivel quanto à probabilidade de ocor- 
rência média anual das várias classes de estabilidade, a 
qual se pode considerar aceitável se os desvios calcu- 
lados numa base anual se situarem numa faixa da ordem 
dos 5— 10%. Admitindo desvios desta ordem de gran- 
deza, o estabelecimento do valor médio anual com base 
em 3-5 anos é geralmente suficiente. Um período desta 
ordem é também o habitualmente considerado nas situa- 
ções de bem maior responsabilidade e implicações que 
são a instalação de Centrais Nucleares. Sem que para 
esta aplicação especifica se considere aceitável a carac- 
terização usualmente feita, nem a margem de incerteza 
admitida, tal valor dá uma indicação indirecta de um 
certo consenso internacional quanto à tipificação das 
probabilidades de ocorrência média anual, No caso ver- 
tente, foi feito o teste ano a ano e verificada a conver- 
gência para um valor médio que seria esperável se o 
número de anos de observação fosse superior. 


Nos capítulos seguintes apresentam-se os resultados 
obtidos. Por eles se constata que em relação a valores 
médios anuais os valores esperáveis não devem afas- 
tar-se de = 5% em relação aos valores apresentados. 
Tal constitui uma validação indirecta para a significância 
estatística do período tratado, quando se procuram valo- 
res médios anuais. 


Cumpre acentuar o carácter anual das médias refe- 
ridas, Os valores trimestrais apresentam maior dispersão, 
a qual naturalmente iria aumentar se fosse considerado 
o mês, a semana, ou o dia, Insiste-se, porém, em que 
uma estatística com base na semana ou até no mês não 
teria grande significado no actual contexto, em «que se 
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visam actividades cuja contribuição individual não é mar- 
cante e em que interessa sobretudo o efeito da sua 
contribuição global para uma poluição de fundo. Numa 
etapa seguinte, seria bem mais importante a obtenção 
da frequência e a caracterização de situações anómalas 
susceptiveis de provocar acidentes, Esta situação deve 
ser considerada no âmbito da elaboração de planos de 
emergência ou de alerta. De igual modo, a consideração 
de grandes contribuições individuais ou de elevada peri- 
gosidade tem de ser abordada na sua perspectiva própria, 
a qual não é a do presente trabalho, muito embora dele 
possa naturalmente decorrer. 


46 — Esclarecidos os limites de validade e anali- 
sadas as implicações da informação base e do seu 
tratamento, os valores apresentados possuem o rigor 
necessário para a finalidade prevista na sua utilização. 


5. PONDERAÇÃO DAS FREQUÊNCIAS 
HORARIAS DE OBSERVAÇÃO 


51 — Sendo a estabilidade da atmosfera fundamen- 
talmente deteminada pelas trocas de energia entre a 
atmosfera e o solo, o simples movimento do Sol deter- 
mina uma evolução ao longo do dia das classes de 
estabilidade. A este efeito local, associado à permuta 
de energia radiante, vem juntar-se os que resultam da 
evolução sinôptica global, WVerificando-se que cerca de 
10% das estações meteorológicas efectua apenas duas 
observações diárias (às 9 e às 18h), caracterizar a 
Dcorrência estatistica apenas com tais observações seria 
introduzir uma distorção que iria ponderar excessivamente 
as classes de estabilidade associadas a menores fluxos 
de radiação incidente, 

Por outro lado, na ausência de método mais rigoroso, 
o qual seria ele próprio incompativel com a precisão 
inerente à tipologia utilizada, a simples consideração das 
escalas globais de evolução sugere que entre os periodos 
considerados a evolução das classes siga evolução seme- 
lhante à que se verifica noutras estações. Pode assim 
estabelecer-se uma regra interpolatória verificando a rela- 
ção existente entre a ocorrência estatística média anual 
que se obteria para uma estação dada se para a mesma 
fossem consideradas 2, 3, 8 ou 24 observações, Tomando 
para ilustração o ano de 1970 e as estações que efec- 
tuam pelo menos B observações diárias e que são: 


QUADRO Il 
575 — Bragança 
532 — Sintra 
536 — Portela 
554 — Faro 
568 — Penhas Douradas 
549 — Coimbra 


571 — Portalegre 
545 — Pedras Rubras 


apresentam-se no quadro seguinte e por classe de esta- 
bilidade, a percentagem de ocorrência anual que se 


32 


obteria para cada uma das estações se em vez das 
8 (ou 24) observações diárias fossem apenas conside- 
radas 3 ou 2. As horas escolhidas para as 3 e as 2 obser- 
vações diárias correspondem ao horário tipico das esta- 
ções que efectuam tal número de observações durante o 


QUADRO Ill 


Valores obtidos para as frequências anuais (1970) 


Estação | C. E. |24obsídia | 8 obs/dia 3 obsídia | 2 obsídia 


95 28 42 


1 
2 12.82 18.58 20.96 
3 13.72 16.04 14.02 
575 4 42.42 44.25 53.40 
5 7.30 6.98 3.26 
6 12.48 8 96 5.95 
7 10.32 4.91 1.98 
1 66 77 1.29 1.94 
2 4.83 5.29 6.55 819 
3 10.33 10.68 14.58 11.39 
545 4 53.41 5117 53.87 67.22 
5 10.64 10.96 7.66 417 
6 15.34 16.35 12.36 5.69 
7 418 | 480 3.69 1.39 
1 95 0.0 
2 943 11.89 
| 3 14.39 14.07 
568 | 4 53.93 63.48 
5 8.48 419 
6 7.81 3.35 
7 3.02 
1 .B4 
2 10.74 
q 12.97 
549 4 64.44 
5 4.88 
6 5.02 
7 É 112 
I AO 14 
2 13.81 13.57 
A 10. 17.52 12.29 
571 4 48.52 43.90 61.29 
5 12.79 1.71 | 5.29 
6 14.57 11.62 | 6.57 
7 2.84 1.33 86 
1 413 419 6.25 
2 9.14 9.85 14.34 
3 21.40 19.50 | 21.32 
532 4 60.20 64.57 55.88 
5 75 42 | 74 
6 350 | 1.05 1.47 
7 .88 42 0.0 
1 36 52 0.0 00 
2 394 378 4.91 4.46 
3 9.50 | 948 14.35 12.40 
536 4 59.89 58.80 66.48 73.54 
5 12.42 io | 333 543 
6 11.75 1215 5.46 3.48 
7 | 44 2.50 46 70 
1 44 3 49 69 
2 6.99 6.46 10.35 13.03 
3 14.59 | 1475 20.15 16.19 
554 4 4091 | 3994 40.84 56.10 
5 811 | 842 8.33 4.66 
6 20.63 21.14 15.66 6.88 
7 


8.34 8.97 421 | 247 
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dia. Apresentam-se também, no mesmo quadro, os valores 
que se verificam no caso das observações serem horá- 
rias, para as estações em que tal informação existe. 


5.2 — Como se pode verificar o aumento do número 
de observações tende a aumentar a percentagem de 
ocorrência das classes intermédias, 2, 3 e 4, que são 
também as mais frequentes, o que é um resultado espe- 
rével. As classes extremas apresentam uma dispersão 
naturalmente maior o que, associado à sua menor per- 
centagem de ocorrência, origina erros relativos maiores. 
Atendendo a que as classes 6 e 7 são classes nocturnas, 
cujas caracteristicas dispersivas originam ampla margem 
de incerteza como anteriormente referido, e à muito baixa 
percentagem de ocorrência da classe 1 (excepto no caso 
de Sintra, devido a caracteristicas localizadas), o esquema 
interpolatório afigura-se como aceitável para as clas- 
ses 2, 3 e 4 Como a classe 5 apenas existe na carac- 
terização de Turner e corresponde a um caso intermédio 
entre as D e E de Pasquill-Gifford e portanto numa zona 
de transição para as classes nocturna, a sua ponderação 
não teria grande significado. No esquema interpolatório 
que finalmente se retêm, para ter em conta as estações 
com apenas duas e três observações diárias, apenas se 
têm em conta as classes 2, 3 e 4, Desta interpolação 
resulta um residuo que se reparte naturalmente pelas 
restantes. A incerteza resultante, devido à baixa frequên- 
cia com que ocorrem, é aceitável em termos da caracte- 
rização global que se pretende. Todavia, deve novamente 
sublinhar-se que no caso de poluidores de elevada peri- 
gosidade ou de grande contribuição individualizada, é 
imprescindivel ter em conta tais situações através de 
estudo local, pois a baixa probabilidade de ocorrência 
anual não significa, em tais casos, ausência de situações 
graves. 


6. VALORES MÉDIOS ANUAIS 
EM GRANDES ZONAS DO TERRITÓRIO 


61 —No anexo | apresentam-se para 1970, 1971, 
1972 e 1973 e cada estação os valores trimetrais e anuais 
de ocorrência das várias classes de estabilidade, 

Na determinação de valores médios anuais a tomar 
como referência tem interesse inferir o limite estatístico 
que se obteria aumentando o número de anos conside- 
rado. Para tal existem vários métodos possíveis, O mais 
simples e intuitivo consiste em verificar as alterações que 
se produzem quando se toma a média de um ano: a 
média da média dos dois, etc. Nos quadros seguintes 
apresentam-se os resultados obtidos para as principais 
estações. Para eles se verifica que os desvios se situam 
nos = 5%. Tal é esperável para as médias anuais e 
justifica a prática internacional corrente de, para este 
tipo de caracterização, se utilizarem sérias de medidas 
com durações semelhantes. 


6.2 — À impossibilidade prática de caracterização 


pontual de todo o território através de medidas directas 
confere grande interesse imediato à determinação de 
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grandes zonas de caracteristicas semelhantes quanto à 
probabilidade de ocorrência das várias classes de esta- 
bilidade com base na informação existente, 

Com os dados existentes e sem recurso a modelos 
fisico-matemáticos muito elaborados e onerosos, a possi- 
bilidade de estabelecer tais zonas reside na admissibili- 
dade de usar um processo interpolatório entre estações. 
A viabilidade de um esquema interpolatório decorre da 
própria natureza do fenômeno em estudo e da localização 
e frequência das observações. 

Tendo em conta as escalas espacio-temporais envol- 
vidas, a natureza da caracterização pretendida, e a loca- 
lização e frequência das observações directas, uma inter- 
polação relativamente simples é adequada muito embora 
se lhe possam levantar objecções decorrentes sobretudo 
do facto de não poder ser directamente validada com a 
informação existente ou de não ser estabelecida com 
recurso a muito maior sofisticação do método, Trata-se, 
porém, de uma solução de compromisso, cujas vantagens 
ultrapassam as eventuais limitações para a finalidade 
prevista. Acentua-se, todavia, que a solução adoptada 
tem carácter provisório e deverá ser reexaminada quando 
for possvel efectuar estudos base mais aprofundados. 
Tais estudos são, aliás, imprescindíveis na consideração 
de actividades com elevado indice de perigosidade, seja 
pela natureza dos efluentes seja pelas quantidades emi- 
tidas, 


6.3 — Com base nos valores médios e tendo em 
conta as observações anteriores, foram traçadas para 
cada ano e classe de estabilidade as isolinhas corres- 
pondentes a frequências iguais de ocorrência O traçado 
das isolinhas utiliza um programa de cálculo automático 
desenvolvido no C.T.A.M.F.U.L. e considera a rede 
interpolatória constituda por triângulos. Tal rede, através 
da sua definição, permite um melhor ajustamento ao peso 
relativo das várias estações, 

Da análise das regiões delimitadas pelas isolinhas 
e da sensibilidade das mesmas a modificação na trian- 
gulação base do ano e da classe de estabilidade resulta 
uma zonagem do território cuja significância se pode 
inferir. 

Nas figuras seguintes apresenta-se, para cada classe, 
a zonagem que resulta ao serem consideradas as médias 
de quatro anos. 

Na representação apresentada, suprimiram-se as iso- 
linha intermédias e as várias zonas apresentam contornos 
propositadamente esbatidos para que do traçado se 
não infira um rigor superior ao fisicamente admissivel. 
Aliás, a simples comparação, ano a ano, da zonagem 
correspondente, mostra (como seria de esperar) altera- 
ções sensíveis. A titulo ilustrativo apresentam-se os 
traçados para a classe neutra correspondente a cada 
um dos anos considerados, Estas variações, que se tor- 
nam imediatamente perceptiveis com a representação 
feita, visualizam simultaneamente a variabilidade espe- 
rável em relação a um ano médio e acentuam a impor- 
tância dos estudos locais pormenorizados sempre que se 
trate de instalações muito importantes. 
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Estação 


575 
Bragança 
(8 obs.) 


1 

2 

543 3 
Viana do Castelo 4 
| (4 obs.) 5 
6 

7 

1 

2 

566 3 

Vila Real 4 
(4 obs.) 5 

6 

7 

1 

2 

545 3 
Pedras Rubras 4 
(24 obs.) 5 

6 

7 

1 

2 

568 3 
Penhas Douradas 4 
(8 obs.) 5 

6 

7 

1 

2 

549 3 
Coimbra 4 

(8 obs.) 5 

6 

7 


540 
Monte Real/B, A. 
(6 obs.) 


40 01 f Co my = 


530 
Cabo Carvoeiro 
(4 nbs.) 
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QUADRO IV 


Frequência de ocorrência das classes de estabilidade 


Estações meteorológicas sinópticas 


1970+ 1971 
MEDIA 


B9 92 
12.75 12.78 
13.43 13.57 
43.99 43.20 

1.44 RT 
12.32 12.40 
918 9.75 
1.05 1.35 
8.11 8.72 
17.20 16.38 
40.84 40.59 
4.76 4.49 
11.75 11.93 
16.52 
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197041971 


Estação C. E MÉDIA 1972 
1 1.24 1.39 
2 7.93 7.34 
571 3 11.45 10.51 
Portalegre 4 48.56 53.28 
(8 obs.) 5 14.27 13.02 
5 13.80 12.06 
7 2.11 240 
1 3.45 1.67 
2 7.97 6.59 
532 3 21.179 16.70 
Sintra/Granja 4 63.35 62.13 
(8 obs.) 5 .69 4.39 
6 2.13 5.80 
7 63 2.i2 
1 5a. 29 
2 SO | 2.11 
536 3 9.89 9,19 
Portela 4 57.58 51.33 
(24 obs.) 5 12.98 | 13.84 
6 12.80 10.08 
7 2592; | 2.56 
1 3.49 2.48 
2 8.22 5.85 
534 3 17.93 17.22 
Montijo/B. A. 4 51.01 58.89 
(5 obs.) 5 888 | 8.22 
6 6.69 5.52 
7 Sis | 1.82 
1 Br | 91 
2 5.15 5.02 
3 8.97 9.59 
4 58.36 62.42 
5 15,24 14.72 
6 8.85 6.60 
7 98 4 

1 197 | 1.31 | 
2 7.64 | 7.13 
554 3 13,11 13.22 
Beja 4 43.86 
(8 obs.) 5 14.63 
5 16.46 
7 2.80 
1 31 
2 5.65 
554 3 13.56 
Faro 4 45.04 
(24 obs.) 5 10.60 
6 18.23 
7 6.61 
1 1.62 
2 3.41 
538 3 9.47 
Sagres 4 69.35 
(8 obs.) 5 8.16 
6 5.03 
7 2.96 
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QUADRO V 


Frequência de ocorrência das classes de estabilidade 
Estações do Continente com 2 e 3 obs./diárias 


1970 + 1971 


Estação GE: MÉDIA 


4.63 
19.09 
23.97 
34.20 

2.42 

5.78 

9.97 


1.03 
10,65 
29.41 
49.11 

4.14 

5.43 

3,64 


AU 
14.24 
15.35 
57.41 

2.39 
7.03 
3.51 

14 
15.80 
19.33 
47.98 


17 
Chaves 
(2 obs.) 


posse 


121 


S. Pedro de Muel 
(2 obs.) 


170 
Setubal 
(2 obs.) 


247 


V. do Alentejo 
(2 obs.) 


“iDunblolrn-|is00bLbr+-Iisonturn-+I40nNbuwurlr- ls: = 


26 


Pedras Salgadas 
(3 obs.) 


114 
Barra do Mondego 
(3 obs.) 


279 
Praia da Rocha 
(3 obs.) 


Jsontorsalsontuwuna lamento 


36 TÉCNICA 460 


E A MAL, E. Lo To AM. FM, E 


distribuição de probabilidade de ocorrência distribuição de probabilidade de ocorrência 
(1970-1973) (1970-1973) 


classe B classe € 
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distribuição de probabilidade de ocorrência 
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Co TAM. FP. Vil, C. TAS MM. Fold. Lo 


distribuição de probabilidade de ocorrência distribuição de probabilidade de ocorrência 


(1970) (1971) 
classe D classe D 
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CT. AM. Pl. L. 


distribuição de probabilidade de ocorrência 
(1972) 


classe D 


E. TLA-M. Fl, L. 


distribuição de probabilidade de ocorrência 
(1973) 


classe D 
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7. CONCLUSÃO 


7.1 —É habitual em estudos deste tipo omitir ou 
diluir limitações que em todas as conclusões existem. 
As entidades responsáveis fomentam usualmente tal ati- 
tude por entenderem que, face à opinião pública, o 
invocar dos resultados de um estudo cientifico ou técnico 
é argumento objectivo de importância decisiva. Reputa-se 
tal atitude como susceptível de gerar um ciclo vicioso 
de graves consequências no presente e sobretudo no 
futuro, Em Portugal, tal atitude tem sido importante fonte 
de inadmissível subserviência a consultores estrangeiros 
em inúmeros dominios. No caso vertente, o sistemático 
assinalar das incertezas e das limitações é não só apli- 
cável ao estudo feito como sobretudo ao modo como 
internacionalmente se tem legislado no capítulo da polui- 
ção atmosférica. Efectivamente, quando se analisa o 
substrato cientifico da legislação existente nos países 
cientifica e tecnicamente evoluídos no concernente a 
chaminés, dispersão de poluentes, etc., a conclusão ine- 
quivoca é o da ausência de adequada fundamentação 
para um número importantissimo de casos onde o mesmo 
se aplica. O facto de a fundamentação ser insuficiente 
não é motivo para que se não introduza regulamentação 
disciplinadora. É motivo sim para que a mesma seja 
continuamente melhorada e sobretudo não ignore os casos 
particulares que a legislação genérica é insusceptivel de 
contemplar. 

7.2 — No âmbito da dispersão de poluentes na atmos- 
fera, houve neste trabalho a preocupação de acentuar 
as limitações intrinsecas da metodologia usada interna- 
cionalmente Tais limitações são aceitáveis tratando-se de 
legislação genérica aplicável a unidades de pequeno 
porte e de baixa perigosidade, E ainda aceitável quando 
a incerteza relativa à quantidade de poluentes emitidos 
não justifica maior rigor. Neste caso, porém, não pode 
iludir-se a gravidade da causa com a incerteza do método, 
Permitir uma incerteza grande nas emissões e devido a 
tal incerteza usar apenas métodos grosseiros na pre- 
visão transforma-se em ciclo vicioso que por vezes se 
torna impossivel de quebrar sem custos sociais muito 
elevados. 

Um dos exemplos mais típicos de tal situação encon- 
tra-se na emissão de produtos radivactivos em Centrais 
Nucleares e na instalação de grandes complexos indus- 
triais fortemente poluidores, como pode ser entre nós o 
caso de Sines e é já o do Barreiro-Seixal. Chama-se 
veementemente a atenção para o contrasenso que será 
a dispendiosa extrapolação de métodos desconhecendo 
as suas limitações, em vez de analisadas estas se 
ultrapassarem as incertezas que introduzem. 

7.3-— As observações anteriormente feitas, devida- 
mente entendidas, permitirão uma fundamentação clara 
da legislação proposta e fixarão os limites da sua razoá- 
vel aplicabilidade. O reconhecimento das incertezas e 
limitações deve ser entendido como estimulante para as 
ultrapassar, não como motivo para que indefinidamente 
se adie a legislação que é urgente promulgar, ou para 
que criticamente se copiem regulamentos estrangeiros 
onde as limitações aqui referidas se desconhecem ou 
foram escamoteadas, 

Deve, finalmente, acentuar-se que para o objectivo 
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inicialmente proposto, a metodologia adoptada é, sem 
dúvida, a mais adequada. Abusivo seria extrapolá-la 
(como é correntemente feito) a situações que a natureza 
do método necessariamente exclui, 

7.4 -— Destinando-se sobretudo a caracterização da 
possibilidade de ocorrência das várias classes de estabi- 
lidade atmosférica à determinação quantitativa das suas 
propriedades dispersas e estas por sua vez ao estabele- 
cimento de fórmulas para o cálculo da dispersão, este 
trabalho é apenas um passo importante nessa direcção. 
Todavia, ele é em si mesmo autónomo pois outras utili- 
zações são possíveis além das referidas. A sua aplicação 
ao estabelecimento de fórmulas de cálculo a utilizar em 
Portugal é por isso objecto de um tratamento separado. 
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NUMERO 460 OUTUBRO 1980 


ANO LV 


VOLUME XLII 


Regime variável em pressão 
Notícia sobre bibliografia recente 


1 — INTRODUÇÃO 


Nos domínios da Hidráulica e da Mecânica dos 
Fluidos, o regime variável em pressão tem merecido, 
nos últimos anos, um renovado interesse por parte de 
investigadores e editores. Este interesse pode ser jus- 
tificado pela ocorrência simultânea das seguintes ra- 
zões: experiência já amadurecida na exploração intensa, 
durante os últimos vinte anos, de modelos computa- 
cionais automáticos de simulação; elevação dos niveis 
de segurança ou de complexidade de operação exigi- 
dos em alguns sistemas elevatórios, nomeadamente os 
circuitos de centrais nucleares e hidroeléctricas com 
bombagem; maior exigência de economia no dimensio- 
namento de sistemas elevatórios de saneamento básico 
a de rega, cuja importância e número têm aumentado 
a nível mundial. Assim, a par de uma tendência natu- 
ral e lógica no sentido da divulgação de novos métodos 
computacionais, o que ocorre periodicamente em todos 
os domínios das ciências aplicadas, verifica-se existir, 
na actualidade, grande procura e interesse nesta ma- 
téria por parte dos profissionais de engenharia, 

É importante referir que nos últimos dez anos têm- 
-sa verificado significativos avanços na modelação ma- 
temática de regimes transitórios em  escoamentos 
difásicos (mistura de fases gasosa e líquida) motivados 
pela ocorrência de pressões subatmosféricas. Nestes 
avanços assumiram especial importância os potentes 
meios de cálculo postos à disposição dos investiga- 
dores. 

A vasta bibliografia sobre o regime variável em 
pressão (ver [1]) foi recentemente enriquecida com 
seis obras importantes, algumas das quais podendo 
vir a constituir marcos de referência tão importantes 
como são actualmente as obras de Allievi (1903), de 
Bergeron (1950), de Parmakian (1955) e de Streeter 
(1967). Nesta nota, apresenta-se uma breve descrição 
dos novos livros com o Intuito de orientar os even- 
tuais leitores interessados na escolha de um texto de 
consulta referente a algum dos tipos de regime variá- 
vel (ver em 6 a terminologia básica). 


Manuscrito recebido para publicação em 5/5/80. 
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ANTÔNIO BETÂMIO DE ALMEIDA 


2 — BIBLIOGRAFIA RECENTE 


A partir da década de 60, regista-se um incremento 
no número de artigos publicados em revistas técnicas 
e de comunicações a congressos, relativos aos regimes 
transitórios. Reuniões especializadas sobre este domi- 
nio passaram a ser realizadas com frequência em paí- 
ses da Europa (incluindo Portugal [2]) e nos E. U. A. 
Em [1] apresenta-se uma síntese seleccionada das pu- 
blicações referentes a este domínio e publicadas até 
1978. Nos últimos anos da década de 70 a bibliografia 
do regime variável em pressão foi acrescida com seis 
importantes livros, em lingua inglesa. Por ordem de 
recepção, estes livros são seguidamente indicados: 


[J-1] Fluid transients in hydro-electric engineering 
practice, de Charles Jaeger, edição de Blackie, 
1977 (413 pág.). 

Hydraulic analysis of unsteady flow in pipe 
networks, de J. À. Fox, edição MacMillan, 1977 
(216 pág.). 

Fluid transients, de E. B. Wylie e V. L. Streeter, 
edição da Mc Graw-Hill, 1978 (384 pág.). 
Applied hydraulic transients, de M. H. Chau- 
dhry, edição de Van Nostrand, 1979 (503 
pág.). 

Modern analysis and control of unsteady flow 
in pipelines, de G. Z. Watters, edição Ann 

Arbour, 1979 (251 pág.). 

Control and suppression of pressure surges in 

pipelines and tunnels, de A, R. D. Thorley e 
K. J. Enever, edição da C. |. R. |. A. 1979 
(112 pãg.). 


LF-2] 


[5-3] 


[C-4] 


LW-5] 


[7-6] 


Seguidamente apresenta-se uma breve descrição 
destas obras as quais serão designadas pelo símbolo 
que antecede os respectivos títulos, 


3 — BREVE DESCRIÇÃO DOS NOVOS LIVROS 
Desde já convém distinguir o livro [J-1] dos res- 


tantes. Com efeito, este livro, constituindo uma com- 
pilação de diversas obras escritas pelo respectivo au- 
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tor, encerra de modo admirável o periodo correspon- 
dente à época anterior ao advento do cálculo automá- 
tico. Ao longo de três partes fundamentais (respeitan- 
tes a oscilações em massa, golpe de ariete e a casos 
concretos de aplicação, respectivamente) e de 38 capí- 
tulos, o livro sintetiza as contribuições da primeira 
metade do século referentes ao comportamento de cha- 
minés de equilíbrio (cálculo de oscilações e da esta- 
bilidade), ao golpe de ariete (métodos de Allisvi e de 
Schnyder-Bergeron), à regulação de turbinas e aos fe- 
nómenos de ressonância e de rotura da veia líquida. 
Ensaios em protótipo e alguns acidentes são descritos 
com pormenor. À utilização do computador e das res- 
pectivas técnicas são apresentadas muito sumariamente 
em dez páginas. O conjunto de referências bibliográfi- 
cas clássicas apresentado é notável, tanto pelo número 
como pela qualidade, realçando-se as referências corres- 
pondentes ao comportamento das chaminés de equili- 
brio em aproveitamentos hidroeléctricos. Com uma 
apresentação luxuosa, o livro [J-1] será no futuro uma 
obra de consulta esporádica mas capaz de substituir 
a consulta de muitas dezenas de referências bibliográ- 
ficas. Nesta nota não se farão mais referências a [J-11. 

Os restantes livros podem considerar-se como re- 
presentativos da experiência adquirida com os métodos 
de cálculo actuais dando os respectivos autores grande 
ênfase aos modelos de cálculo automático baseados 
no método das características. Em [F-2] e [5-3] en- 
contram-se ainda breves referências ao método gráfico 
mas em [C-4] e [W-5] toda a formulação é orientada 
para o uso do computador digital. 


Todos estes livros apresentam as bases teóricas 
fundamentais do regime variável em pressão e do mé- 
todo das características bem como indicações sobre a 
modelação matemática de turbomáquinas hidráulicas e 
dae alguns dispositivos de protecção. São também apre- 
sentadas introduções aos regimes variáveis em escoa- 
mentos com superfície livre ([F-2], [5-3] e [C-4]), 
em fluídos gasosos ([F-2] e [5-3]) e em outros do- 
mínios especiais ([5-3]). No que concerne os sistemas 
hidroeléctricos os livros [5-3] e, em especial, [C-4] 
apresentam capítulos de interesse referentes a: regu- 
lação de turbinas, ressonância ([5-3] [C-4]) e cha- 
minés de equilíbrio ([C-4]). Em [C-4] existem capí- 
tulos especiais respeitantes a regimes transitórios em 
circuitos de centrais nucleares e em oleodutos. 

Os livros [F-2], [W-5] e [7-6] são orientados 
para a análise do regime variável em sistemas eleva- 
tórios equipados com turbobombas. As bombas de êm- 
bolo e a determinação de leis de fechamento adequa- 
das para válvulas de controlo em condutas gravíticas 
merecem capítulos especiais em [5-3]. Exceptuando 
o caso de [F-2], todos os restantes livros apresentam 
programas de cálculo automático para a análise de 
sistemas simples. 

Não obstante todos os livros apresentarem infor- 
mações gerais sobre dispositivos de protecção, é em 
[7-6] que o projectista pode encontrar um conjunto de 
elementos úteis para o dimensionamento de alguns 
desses dispositivos. 


No que respeita a cavitação e a rotura da veia 
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líquida, [F-2], [S-3] e [C-4] dão indicações e refe- 
rências bibliográficas relativas a avanços recentes na 
modelação matemática de escoamentos difásicos. 


4 — CONSIDERAÇÕES FINAIS 


Uma apreciação qualitativa dos livros referidos 
nesta Nota terá de ter em consideração os diferentes 
tipos de leitores que recorrerão aos mesmos. Ássim, 
o analista de sistemas hidráulicos em pressão já fami- 
liarizado com modelos de cálculo automático, poderá 
encontrar em três das seis obras ([F-2], [5-3] e 
[C-4]) algumas sugestões e pistas novas para a reso- 
lução de casos complexos. Para quem desejar uma 
iniciação muito rápida na modelação de sistemas ele- 
vatórios e estiver pouco preocupado com considerações 
teóricas, encontrará em [F-2] e [W-5] bons elementos 
de referência. Finalmente, o projectista de condutas 
que dispense pouco tempo no aperfeiçoamento de téc- 
nicas de análise irá encontrar em [T-6] algumas res- 
postas positivas. Não se quer deixar de referir que, 
ao contrário do que sucedeu com os métodos de cál- 
culo anteriores, o domínio das técnicas recentes de 
computação não pode ser conseguido pela simples lei- 
tura de um livro. Com efeito, entre os esquemas e 
algoritmos de cálculo apresentados e os programas de 
cálculo operacionais, para simulação de casos reais, 
existe um grande fosso que só a experiência adquirida 
numa prática constante de modelação e no acompa- 
nhamento das testagens dos respectivos programas po- 
derá vencer. 

Todas as obras apresentadas têm, para os técnicos 
portugueses, o inconveniente de serem escritas em in- 
glês. Anuncia-se, para breve, contudo, o aparecimento 
de algumas publicações em lingua portuguesa sobre 
este domínio. 


5 — REFERÊNCIAS BIBLIOGRÁFICAS 


[1] Almeida, A. B.— O golpe de ariete em condutas 
elevatórias. Bibliografia e sintese dos conhecimen- 
tos. LNEC, 1979. 

[2] Almeida, A. B.--— O golpe de ariete em condutas. 
Textos apresentados no Seminário com o mesmo 
título. LNEC, 1979, 


6 — TERMINOLOGIA BÁSICA 


OU regime variável em pressão («unsteady pressure 
flow») é aquele em que a pressão e a velocidade do 
escoamento variam com o tempo e ao longo das con- 
dutas. Nestas variações, que apresentam características 
ondulatórias («pressure surgesr), não são usualmente 
incluídas as flutuações locais motivadas pela turbulên- 
cia. Designam-se por regimes transitórios («fluid tran- 
sientsy) os regimes variáveis que ocorrem na transição 
entre dois regimes permanentes. No caso de um regime 
transitório mobilizar as forças elásticas do líquido e 
das condutas, é usual designá-lo pelos termos suges- 
tivos de golpe de ariete ou choque hidráulico (coup de 
bélier ou waterhammet). 
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Notas sobre a teoria da semelhança 


RESUMO 


Após discussão do problema da invariância das 
equações da Física nas mudanças de sistemas de unida- 
des é apresentada uma definição de semelhança física. 

É em seguida apresentado o problema das possibili- 
dades de realização de sistemas físicos semelhantes. 

Uma notação proposta neste texto permite a reso- 
lução quase automática dos problemas de semelhança 
(das relações entre as grandezas físicas de sistemas 
semelhantes). 

Estas notas foram redigidas para serem usadas pelos 
alunos da cadeira de Introdução à Fisica dos Meios 
Continuos do |, 5, T. Com a sua publicação na «Técnica» 
pretende-se contribuir para que nas cadeiras de Fisica 
das escolas de Engenharia passe a ser dada maior aten- 
ção a este assunto e a este tipo de problemas, 


APRESENTAÇÃO 


A diversidade dos sistemas de unidades em uso 
até há bem pouco tempo obrigou a que na formação 
dos técnicos fosse olhado com particular cuidado e 
atenção o problema da alteração dos valores numéricos 
das grandezas físicas nas mudanças de sistemas de uni- 
dades. 

Gerações de estudantes de Fisica e de Engenharia 
(e desde há algumas décadas do secundário) foram 
assim, por via do problema da mudança dos sistemas 
de unidades, sensibilizados para o problema da dimen- 
são, da homogeneidade e da invariância das fórmulas 
fisicas. 

As duas questões, porém, não se confudem. A pri- 
meira insere-se na segunda sem a esgotar. 

Saber mudar de sistema de unidades é, com efeito, 
um mero problema técnico (quase vazio de conteúdo 
ísico) facilmente resolvido por quem domine minima- 
mente o problema da dimensão das grandezas físicas. 

O segundo problema, o da homogeneidade e dimen- 
são das fórmulas físicas e do significado profundo da 
sua invariância, é, pelo contrário, uma questão muito 
delicada, que se relaciona com a construção global da 
Física e da própria Geometria (encarada aqui como o 
primeiro capítulo da Fisica). Não se trata, de modo al- 
gum, de uma questão adicional que surge em determi- 
nada altura, mas sim de algo que está implícito nas 
nossas primeiras noções. 


Manuscrito recebido para publicação em 20/10/80. 
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ABSTRACT 


After discution of the problem of the invariance of 
Physical equations in the changes systems of unit, a 
difinition of physical similarity is presented. 

Then the problem of possibility of realization of simi- 
lar physical systems is presented. 

A notation proposed in this text gives the possibility 
of an almost automatic resolution of the problems of 
similarity (of the relations between physical variables of 
similar systems). 

These notes were written to be used by scholars 
of «Introduction to the Physics of Continous Media» at 
the |. S. T. With the notes here presented we are trying 
to given a contribution so that in several engineering 
schools a greater attention to this subject and this kind 
of problems shall be given. 


Adicionalmente trata-se de um assunto relacionado 
com um problema de grande importância técnica, em 
especial no campo da engenharia, o problema do com- 
portamento e das possibilidades de reprodução de si- 
tuações semelhantes em escalas diferentes. 

Os problemas de semelhança (em que há que con- 
siderar sistemas físicos diferentes) são algo mais com- 
plexos do que Os simples problemas de unidades (em 
que temos um só sistema físico e nos limitamos a mudar 
de sistema de unidades). 


Infelizmente, no ensino tradicional, os segundos 
têm ofuscado os primeiros. 

Em consequência, é corrente ver alunos que doml- 
nam bem o problema das unidades ficarem desorientados 
quando se lhes passa um problema do tipo: 


«Uma mola vibra com uma determinada frequência f,. 
Com que frequência vibra uma mola semelhante, feita 
com o mesmo material, e de dimensões lineares 2 vezes 
maiores?» 

Estando a diminuir a importância prática dos pro- 
blemas de mudanças de unidades por ser cada vez 
menor o número de sistemas de unidades em uso, 
parece ter chegado a hora de os substituir nos cursos 
de Física por problemas de semelhança. E não nos ter 
chegado ainda esta tendência trazida por livros estran- 
geiros não é razão para ficarmos à espera. 
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As notas que se seguem são uma versão corrigida 
de um texto em uso na cadeira de Introdução à Fisica 
dos Meios Contínuos do |, S. T., desde há dois anos, 

Por Teoria da Semelhança entendem-se aqui todos 
os assuntos relacionados com os problemas de seme- 
lhança, o que, necessariamente, engloba as questões 
incluídas no atraz designado segundo problema. 

A Teoria da Semelhança não é, como dissemos, um 
capítulo autónomo, susceptivel de ser ensinado como 
matéria adicional em relação a capitulos anteriores de 
Fisica. 

Sendo assim, a melhor maneira de levar um aluno 
a dominar o assunto é fazê-lo reflectir sobre o modo 
como aparecem os primeiros resultados e se escrevem 
as primeiras fórmulas da Fisica (e da Geometria). 

Estas notas, organizadas como um encadeado de 
comentários, definições e problemas, destinam-se, fun- 
damentalmente, a estimular essa reflexão. 


Espera-se que os leitores estejam familiarizados com 
as dimensões das grandezas físicas e fórmulas de trans- 
formação dos seus valores nas mudanças de sistemas 
de unidades, matéria hoje, dum modo geral, assimilada 
no decorrer dos estudos secundários. 


A esses leitores se pretende mostrar que para 
resolver um problema de semelhança (no âmbito de um 
capítulo da Fisica conhecido) não é precisa nenhuma 
matéria adicional mas somente uma reflexão sobre o 
que já se sabe, 

O uso de uma notação adequada permite-nos apre- 
sentar um processo quase automático de resolução dos 
problemas. Pretende-se, naturalmente, que esse automa- 
tismo seja um automatismo compreendido. ; 


Estas notas atingirão plenamente o seu objectivo 
se ajudarem os leitores a ver, face uma situação física, 
as possíveis situações físicas semelhantes, 


|— COMENTÁRIOS GERAIS 


1 — Invariância das fórmulas e semelhança de sistemas 
físicos 


As leis da Física — olhadas como peças básicas da 
nossa compreensão do Universo — são necessariamente 
aplicáveis a escalas diferentes (a intervalos de variação 
mais ou menos largos das grandezas em Jogo) sem o 
que seriam simples registos de resultados particulares. 
Por isso se exprimem por expressões matemáticas em 
que as grandezas físicas figuram como variáveis a que 
correspondem, em cada caso concreto e uma vez esco- 
lhidas as unidades, determinados valores numéricos. 


A experiência revela que, pelo menos em primeira 
aproximação, vários resultados da Física são indepen- 
dentes da escala (dos valores grandes ou pequenos) 
das grandezas em jogo. 

A atitude primeira do físico é assim a de escrever 
fórmulas não só aplicáveis a diferentes escalas mas 
que traduzam também resultados independentes das es- 
calas, 
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Para satisfazer esta condição as fórmulas devem 
ser homogéneas. 

É o que se passa com as fórmulas Iniciais da Física 
a saber, as fórmulas da Geometria, as fórmulas de 
definição das grandezas cinemáticas, e a fórmula funda- 
mental f=m.a. 

Os sistemas de unidades que adoptamos (depois 
de escolher as grandezas fundamentais e de calcular 
a: dimensões das grandezas derivadas) são construídos 
de modo a assegurar a condição fundamental de todas 
as fórmulas citadas se manterem válidas (terem a mesma 
expressão matemática e neste sentido serem invariantes) 
em todos os sistemas de unidades. 


Do ponto de vista puramente numérico diferentes 
sistemas físicos estudados num mesmo sistema de uni- 
dades podem confundir-se (ser equivalentes dizem al- 
guns) com um único sistema físico estudado em dife- 
rentes sistemas de unidades. 

Esta equivalência, quando se verifica, traduz a possi- 
bilidade de reproduzir a escalas diferentes o mesmo 
sistema físico que, naturalmente, não será o mesmo mas 
sim diferentes sistemas que diremos semelhantes. 

Podemos assim dar a seguinte definição de seme- 
lhança: 

Dois sistemas físicos são semelhantes se, para as 
grandezas que nos interessa considerar, do ponto de 
vista numérico, tudo se passar entre os dois sistemas 
(estudados num mesmo sistema de unidades) como se 
se tratasse de um único sistema físico estudado em 
dois sistemas de unidades diferentes (ver notação 
adiante). 

Partindo de um sistema físico (e limitado o nosso 
interesse a um conjunto limitado de grandezas) a toda 
a mudança de sistema de unidades (e podemos Imagi- 
narmos quantas quizermos) corresponde, assim, o con- 
junto dos valores numéricos (no sistema de unidades 
inicial) das grandezas de um sistema físico semelhante 
ao primeiro... no caso deste sistema poder existir, o que, 
como veremos, nem sempre se verifica. 


2 — Grandezas variáveis e constantes universais 


Nas fórmulas iniciais acima referidas figuram unica- 
mente grandezas (variáveis) e números puros (constan- 
tes). No caso de estarem unicamente em jogo fórmulas 
deste tipo é sempre possivel (teoricamente entenda-se) 
dado um sistema inicial, construir todos os sistemas 
semelhantes correspondentes a todas as (infinitas) mu- 
danças de sistemas de unidades que queiramos imaginar. 

Mas tal não é possível em situações mais comple- 
xas quando nas fórmulas figuram as chamadas cons- 
tantes universais. 

Velamos num exemplo simples como aparecem estas 
constantes. 

Admitamos que tinhamos já definido processos de 
medida das grandezas distância, massa e força, que 
tinhamos adoptado um sistema de unidades, o sistema 
c.g.s.. por exemplo, e que começávamos por via expe- 
rimental a estudar a lei da atracção gravitacional. 
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Ao fim de demorados estudos encontrariamos o 
resultado 


mm 
f = 6,673.10%8 


Tinhamos encontrado uma nova lei fisica expressa 
pela expressão anterior no sistema c.g.s. Se quizéssemos 
passar para um outro sistema, para o sistema internacio- 
nal por exemplo, o nosso problema seria um simples 
problema de matemática e facilmente obteriamos 


mm 
f = 6,673. 101 


Apresentadas as coisas deste modo podemos dizer 
que a lei da gravitação se exprime por funções diferen- 
tes das massas e das distâncias nos dois sistemas de 
unidades. Podemos unificar as duas expressões escre- 
vendo 


mm 


e dizendo que G é uma constante universal de dimen- 
são M! Lº T? que toma os valores indicados nos dois 
sistemas. 

Em resumo temos o seguinte: trabalhando num certo 
sistema de unidades encontramos para a lei da gravita- 
ção uma expressão do tipo f=f, (mm',r). Conhecida a 
função f, o problema da passagem para a expressão 
f=f, (m,m',r) válida num segundo sistema de unidades 
é um mero problema de matemática. As funções f, e f, 
são diferentes. O que há de estranho na lei de gravitação 
de Newton é o facto de podermos unificar as duas 
expressões por melo da introdução de uma única cons- 
tante universal. De conseguirmos com uma só constante 
formular uma lei física invariante na mudança de sistemas 
de unidades. 

Note-se que, no estudo de um fenómeno novo, a 
unificação podia exigir a priori a introdução de várias 
constantes universais, A estranha situação encontrada, 
de ser necessária uma única constante, é, porém, a regra 
geral em Fisica em todos os fenómenos fundamentais. 

Limitamo-nos a registar o facto que, naturalmente, 
é um tema de meditação. 


3 — Limitações na construção de sistemas semelhantes 


Encontramos um segundo tipo de leis em cujas 
expressões matemáticas figuram grandezas variáveis e 
constantes universais (para as imaginarmos com valores 
diferentes temos de nos imaginar noutro Universo). 


Vejamos como se põe o problema da semelhança 
nesta nova situação. 


Suponhamos que temos um sistema físico cujo es- 
tudo envolve uma lei física onde figura uma constante 
universal G. Num sistema de unidades inicial as dife- 
rentes grandezas do sistema e a constante G tem deter- 
minados valores. Podemos imaginar livremente diferentes 
sistemas de unidades. Em geral, nesses novos sistemas 
a constante G tem valores diferentes. Não podemos por 
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isso, em geral, imaginar sistemas reais semelhantes ao 
sistema inicial, que correspondam no sentido atrás des- 
crito a todas as mudanças de sistemas de unidades, 
dado que em muitos desses sistemas a constantes G 
deveria ter um valor diferente do valor que tem no nosso 
Universo. Só são pois possíveis os sistemas semelhantes 
que correspondem às mudanças de sistemas de unidades 
em que o valor de G não se altera, 


Consideremos agora um sistema físico cujo estudo 
exija a consideração de um certo número de leis em que 
figurem as constantes universais G,, Ga, ... Gp. Em situa- 
ções deste tipo poderá não ser possível construlr a uma 
escala diferente um sistema físico semelhante ao sistema 
inicial, por ser impossível encontrar dois sistemas de 
unidades distintos em que os valores numéricos das 
constantes universais sejam os mesmos. 

Consideremos agora a situação mais restrita de 
estarmos a estudar um problema particular, de mecã- 
nica dos fluidos por exemplo, e de desejarmos construir 
um sistema semelhante ao sistema estudado mas sem 
mudar de fluido. No caso da densidade do fluido ter 
influência no fenómeno em estudo vemos imediatamente 
que p desempenha neste problema restrito um papel 
idêntico ao das constantes universais nos problemas 
anteriores. Só nos será possivel construir sistemas seme- 
lhantes que correspondam a mudanças de sistemas de 
unidades que não alterem o valor numérico de p. 

Se além de p tiverem influência no problema em 
estudo toda uma série de outras grandezas caracteris- 
ticas do fluido podemos (será a situação corrente) estar 
impedidos de construir um sistema semelhante a uma 
escala diferente com o mesmo fluido. 


Vemos que os problemas de semelhança são mais 
complexos que os problemas de unidades. 

Há um ponto importante a sublinhar. Sendo, em 
princípio, todas as situações físicas influenciadas por 
todas as leis da Fisica nunca é possível reproduzir seme- 
lhantemente, em rigor absoluto, uma situação física a 
uma escala diferente. 


Nos problemas correntes, porém (e se não fosse 
assim não resolveriamos um único problema) estudamos 
situações em que as grandezas em jogo estão relaciona- 
das por um número limitado de leis fisicas (por serem 
desprezíveis todos os outros efeitos). Nestes casos é, 
às vezes, possivel, mas nem sempre como vimos, repro- 
duzir semelhantemente as situações (os sistemas) físicos 
em estudo em escalas diferentes. 


Compete, naturalmente, ao físico ou engenheiro que 
estuda determinada situação conhecer, ou pelo menos 
delimitar, as grandezas e leis físicas que podem influir 
no seu evoluir. É a partir desses elementos que pode 
saber se é possivel ou não reproduzir situações seme- 
lhantes. 


Problemas simples são aqueles em que os fenó- 
menos em estudo obedecem a leis rigorosas em que 
figuram poucas constantes, universais e caracteristicas 
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dos materiais. São deste tipo os problemas que a seguir 
se encontram neste texto. À nosso ver tem o seu lugar 
nas cadeiras de Fisica, 

Problemas mais complicados são aqueles em que só 
se consegue uma semelhança mais ou menos aproximada 
conforme é maior ou menor a influência das grandezas 
cujos efeitos somos obrigados a ignorar. São problemas 
deste género que se põem aos engenheiros que traba- 
lham em modelos reduzidos, 

Problemas extremamente complicados são aqueles 
em que, dada a natureza própria do problema, é neces- 
sário trabalhar com modelos reduzidos, mas estes mo- 
delos se afastam substancialmente das condições de 
semelhança por a multiplicidade dos fenômenos em jogo 
impedir, mesmo eproximadamente, reproduções seme- 
lhantes a escalas diferentes. É o que se passa, por 
exemplo, no estudo dos incêndios em modelos redu- 
zidos. 


ll — EXEMPLOS E REGRAS PRÁTICAS 
1 — Notações 


Usaremos quase exclusivamente sistemas de unida- 
des em que as grandezas fundamentais são o compri- 
mento, a massa, o tempo e a intensidade da corrente 
eléctrica. Representaremos por L, M, Te | as dimensões 
destas grandezas. 


Dada uma grandeza fisica q explicitaremos as suas 
dimensões com a notação tradicional: 


lg] =MºLº TÍ N 


Sendo gq, e g; os valores em dois sistemas de 
unidades SU, e SU, da grandeza q associada a um 
dado sistema físico temos a relação: 


q By 0 
go = 941 Mjz Lj Tiz ly 


em que |, por exemplo, é a relação entre a unidade 
de intensidade do sistema 1 e a unidade de intensidade 
do sistema 2: 


unidade de intensidade de SU, 


2 “unidade de intensidade de SU,| 


92 
Usamos ainda a notação 9, = — 


Vamos introduzir, agora, uma notação adicional ade- 
quada ao estudo dos problemas de semelhança. 


Admitamos que foi possivel estabelecer uma corres- 
pondência entre as grandezas de dois sistemas físicos 
diferentes SF, e SF, (que em particular associe a todo 
o ponto e a todo o instante do primeiro sistema um 
ponto e um instante do segundo). O símbolo 961) 


representará o valor numérico de uma grandeza de SF, 
num dado sistema de unidade, e é PEN valor numérico 


A8 


da grandeza correspondente de SF; no mesmo sistema 
de unidades. 


3 (2) 


(m 3(1) 


Adoptamos ainda a notação g 


FIG. 1 


No caso de termos dois sistemas físicos e dois 
sistemas de unidades usamos as notações 9,1); € 912 
para representar os valores numéricos de uma grandeza 
de SF,, respectivamente nos sistemas de unidades SU, 
e SU, e 9») € 92; Para representar os valores 
numéricos da grandeza correspondente de SF; nos dois 
sistemas de unidades. 


Temos portanto 


9 (2)1 9 (2 
1 
ERR 9 (ni 


9 (192) & 


Em conformidade com a definição anterior de seme- 
lhança, dois sistemas físicos SF, e SF; são semelhantes, 
se for possivel encontrar dois sistemas de unidades SU; 
e SU; tais que, para todas as grandezas que nos inte- 
ressa considerar, se verifique 9,2) = 9(1)2º 


Sendo SF, e SF, dois sistemas semelhantes, entre 
os valores de uma grandeza de dimensão (&, B, Y, 0) 
de SF, e o valor da grandeza correpondente de SF; 
verifica-se a relação 


9 (21 I (nz 


(mi (mi 


AR E 
= My Ly To lo = 


9(12)= 


No caso de qg ser uma distância d usaremos a 
n = = E 
otação 9 (19) Fanta (io) No caso de g ser uma 
massa, um intervalo de tempo, ou uma intensidade de 
corrente, usaremos para g a»! respectivamente, as no- 


tações Ma)! Te» (7 
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Escrevemos portanto 


g 
(2) CG BP 4t.8 
9 (12) Rd M (2) L(42) LR 2) 
Ja) 


fórmula fácil de fixar porque muito semelhante à ante- 
rior, mas em que já não há referência aos sistemas de 
unidades. 


Lao»! Mao! Tao * ta 


zem a relação de semelhança fisica entre os dois sis- 
temas semelhantes SF, e SF.. Assim, se o segundo 
tiver dimensão lineares duas vezes superiores às do 
primeiro temos L(,») = 2 e se o segundo sistema for 
duas vezes mais rápido do que o primeiro temos 
Taz = 0,5 


são números que tradu- 


2 — Caracteristicas dos sistemas semelhantes e sua 
realização prática 


Na primeira parte deste texto preocupamo-nos com 
o problema da possibilidade de existência teórica de 
situações (sistemas) semelhantes. 

Vamos agora debruçarmo-nos sobre o problema prá- 
tico da realização efectiva de sistemas semelhantes a 
um dado sistema que desejamos estudar. É o problema 
que se põe aos engenheiros que trabalham com modelos 
reduzidos. 


Da definição anterior decorre que dois modelos 
semelhantes gozam das seguintes propriedades: 


|! — São geometricamente semelhantes. 


2 — Para cada grandeza (considerada) há proporcio- 
nalidade entre os valores associados aos elementos 
correspondentes (ou pontos homólogos) de um e outro 
sistema, 


3— Os coeficientes de proporcionalidade 9;(12) en- 
tre os valores num e noutro sistema de todas as gran- 
dezas consideradas verificam as relações 


0, 
i(2) + 
ço | Ni | Bi Ti [Oi 


9;(12) & 02» “ax “2 “(2 


em que os expoentes «;, 8, Y; e 8, indicam as dimen- 
sões das referidas grandezas. 

Como podemos assegurar que estas considerações 
se verificam num modelo reduzido? 

Na prática o problema põe-se do seguinte modo: 

Se soubermos que o evoluir do sistema em estudo 
depende de um certo número de leis (cujas equações 
podem, eventualmente, ser desconhecidas, mas cujas 
constantes temos de conhecer) e se à escala do modelo 
o evoluir do sistema depender das mesmas leis; se, 
adicionalmente, certos dados determinarem o evoluir do 
sistema em estudo e for possivel reproduzir no modelo 
reduzido dados semelhantes respeitando as condições 3, 
então, todo o sistema em estudo é reproduzido seme- 
lhantemente no modelo, 
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As constantes (com os mesmos valores, ou com 
alterações fixas, na passagem do sistema para o modelo) 
impõem limitações na variação dos coeficientes M 


Ley Taz q Lazy 


com dificuldades de natureza prática, que condicionam 
a possibilidades de realização de modelos. 


(12)! 
São estas limitações, conjuntamente 


3 — Exemplos 


Consideremos dois sistemas hidráulicos semelhan- 
tes. 

Admitamos que o evoluir dos sistemas é influen- 
ciado pelas constantes: densidade p, coeficiente de vis- 
cosidade L([H]=M L-! T-') e aceleração da gravidade q. 

Se em ambos os sistemas o liquido for a água e se 
os dois estiverem no mesmo planeta as relações: 


= =] -1 RR 
Eao=Meazy lo a =] 
(12) = Ma Liz = 1 

E E = 
9 (12) & L(42) Ta = 


implicam: Log=1;Mag=4Tay= bh vagsl 


Para serem semelhantes os dois sistemas têm de 
ter iguais escalas de comprimentos e de velocidades o 
que significa serem iguais, 


o is 


Admitamos que o primeiro sistema SF, esta na Terra 
e que queriamos construir um sistema semelhante SF» 
com o mesmo líquido, mas num outro planeta com a 
aceleração da gravidade 8 vezes inferior à da Terra. 
Neste caso temos: 


92) = ea) Trio = 8” 

Hay = Meo e Tag =1 € Pay = Ma) Lt») ha 
Estas relações dão-nos: 

bag =? + Mg=8 à Tg=4 ti Voy = 


o que significa que, no outro planeta, o segundo sistema 
deve ser construido numa escala duas vezes maior mas 
que nele as velocidades devem ser duas vezes mais 
pequenas. 

Temos ainda neste caso: 


f M L Te = 4 


(12) (2) 2) (12) 


o que significa que nos dois sistemas as forças (sobre 
elementos correspondentes) são iguais. 


OO 


Admitamos que em vez de mudar de planeta mudá- 
vamos de liquido e iamos usar em SF; um óleo de 
densidade 0,8 e de viscosidade 4 vezes superior à da 
água. 
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Teriamos neste caso: 


8 
Va = 4 E 6 , 8) =| , 


relações que nos permitem chegar a: 


— E1/3 = 90 
9 =º É M (42) 20 


20. 


7 
Lo=5/3=292 , T 


— Bly3 Er 
= 5/3 q 2) E 


“(12) 
Neste terceiro caso as dimensões lineares, as veloci- 
dades e as forças seriam maiores no sistema 2. 


4 — Semelhança de Heynolds 


Vimos que um sistema cujo estudo envolva a con- 
sideração das grandezas d, P, |, g não pode (no mesmo 
planeta e com o mesmo líquido) ter uma reprodução 
semelhante a uma escala diferente. 

Admitamos, porém, que no problema que nos in- 
teressa estudar a influência de g é nula ou diminuta, 
É o que se passa, por exemplo, no estudo do desloca- 
mento de um sólido no interior de uma grande massa 
liquida (ou mesmo de um gás desde que o efeito da 
compressibilidade seja diminuto). 

Neste caso as únicas grandezas que nos interessa 
considerar com valores impostos são p e 


As condições: 


[1 T:! 


Eua gia qa Ma + 


cão o ap 
1 e Pay = so “(toy T 


Hei) = M 


dão-nos a liberdade de adoptar, por exemplo: 


d (3) 


=——— = 10"! 
e dm , 


L 


O que nos conduz a: 


= os My = = 
=102 ; v 10 4 fo =! 


Vemos que para ter um sistema semelhante com 


uma escala 10 vezes inferior temos de ter uma veloci- 
dade 10 vezes superior, o que na prática dificulta 
bastante a realização de sistemas semelhantes deste 
tipo. 


LI] had 


(Em sistemas geometricamente semelhantes, esco- 
lhida convenientemente uma dimensão linear e uma 
velocidade v, é usual definir um número sem dimensões 


designado por número de Reynolds. 

Dois sistemas semelhantes (no sentido da definição 
que atrás apresentamos) têm números de Reynolds 
iguais. No caso de as únicas grandezas em jogo serem 
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p, v, d, w esta igualdade é condição suficiente de seme- 
lhança. 

Este tipo de semelhança (semelhança de Reynolds) 
é como vimos dificil de reproduzir num modelo. Por esta 
razão, nos estudos em modelo reduzido em hidráulica 
e em aerodinâmica raramente se exige (porque seria 
difícil) que no modelo o número de Reynalds seja igual 
ac do original (que haja semelhança entre um e outro 
sistema). Para poder utilizar os resultados obtidos no 
modelo, torna-se por isso necessário um conhecimento 
do modo de como as grandezas (as forças, por exem- 
plo) dependem do número de Reynolds, 


5 — Semelhança de Froude 


Consideremos agora um problema de mecânica dos 
fluidos em que a influência da viscosidade seja despre- 
zivel. Admitamos que usando o mesmo liquido e sem 
mudar de planeta queriamos fazer uma reprodução seme- 
lhante a uma escala diferente. 

Neste caso as duas únicas relações: 


Es -3 = 
Pça) = Mao btt E] 
E É a 
9 (19) L,19) Tão | 


permitem-nos escolher, por exemplo. Lç, = 10! o que 
o que implica: 


-1/2 1/2 


=10 , E. = 103, 


Tal significa que, se no modelo as dimensões linea- 
res forem dez vezes menores, as velocidades devem ser 
raiz de dez vezes menores, as forças mil vezes menores 
e as pressões dez vezes menores, 

Nestes casos (em que a influência da viscosidade 
é nula ou diminuta) podemos com facilidade construir 
modelos semelhantes a escalas reduzidas. 

É o que se fez, por exemplo, no estudo do molhe 
dao Sines. À influência da viscosidade da água na acção 
do mar sobre um molhe é praticamente nula. A gran- 
deza | não figura, com efeito, nas equações que per- 
mitem com uma aproximação razoável estudar o pro- 
blema. As grandezas que influenciam os cálculos são 
unicamente pg, ve d. 

Num modelo é possivel testar a possibilidade do 
mar arrastar blocos sólidos. Com efeito, o equilibrio de 
um bloco depende, fundamentalmente, do seu peso, que 
diminui com o cubo da dimensão linear, e das forças 
hidráulicas que como acabamos de mostrar se reduzem 
do mesmo modo. 

Não é porém fácil testar a possibilidade do mar 
partir os blocos, Com efeito, no modelo reduzido as 
pressões (e portanto também as tensões internas) apa- 
recem reduzidas na mesma escala das distâncias. 


Em Sines o mar partiu os dolos. Para num modelo 
na escala de 1 para 10 fazer ensaios de ruptura dos 
dolos teria sido necessário construir os modelos redu- 
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zidos de dolos num material com a mesma densidade 
dos dolos reais, mas com uma carga limite de ruptura 
1) vezes menor. 


E não seriam ainda considerados outros efeitos 
como os da compressibilidade da água. 


É usual nos problemas de mecânica dos fluidos 
definir um número sem dimensões 


V 


F=—-— 
vg d 

em que v e d são uma velocidade e uma distância 
convenientemente escolhidas (de igual modo em todos 
os sistemas semelhantes) designado por número de 
Froude. 

Em dois sistemas semelhantes o número de Froude 
(tal como o de Reynolds) é o mesmo, A igualdade do 
numero de Froude é, assim, uma condição necessária 
de semelhança. Mas esta condição é também uma con- 
dição suficiente nos problemas em que as únicas gran- 
dezas com importância significativa são, Pp, g ve d. 
Quando é este o caso, e quando os números de Froude 
de dois sistemas são iguais, falamos em semelhança 
de Froude, tal como no caso anterior falávamos em 
semelhança de Reynolds). 


6 — Escolha de tipo de semelhança 


Dado um sistema físico cujo estudo envolva a con- 
sideração de um número elevado de variáveis, vemos 
que, em geral, só é possível construir sistemas seme- 
lhantes com iguais dimensões lineares. 


A construção de modelos semelhantes com escalas 
diferentes só é possivel quando os condicionalismos 
relativos às grandeza que determinam a evolução (ou 


o estado) de um sistema não impõem a relação Lay= 1. 


O estudo dos problemas de semelhança começa, 
assim, pela selecção de quais estas grandezas. Este 
conhecimento não é fornecido propriamente pela Teoria 
de Semelhança, mas sim pelo capitulo ou capitulos da 
Física envolvidos no estudo do sistema em causa. 


São assim, por exemplo, conhecimentos de Mecá- 
nica dos fluidos que nos permitem dizer, face a um 
determinado problema de movimento de fluidos, se a 
reprodução semelhante a escala diferente do sistema 
em estudo tem que ver com a semelhança de Froude, 
a semelhança de Reynolds, ou outros tipos de seme- 
lhança relacionados com outras grandezas determinantes 
na evolução do sistema, Tal, exige, nalguns casos, conhe- 
cimentos bem assimilados. 


Noutros campos da Fisica, porém, a indicação das 
grandezas intervenientes nos fenómenos é elementar. Os 
problemas de semelhança aparecem, assim, como pro- 
blemas simples, de solução quase imediata, que podem 
ser resolvidos com uma metodologia que quase se trans- 
forma num automatismo. 


Todos os problemas da série que a seguir se apre- 
senta são deste tipo. 
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7 — Problemas 


P. 1 — Uma barra encastrada vibra com uma frequên- 
cia de 2 vibrações por segundo. Qual deverá ser a 
frequência de vibração de uma barra semelhante feita 
com o mesmo material, igualmente encastrada, e de 
dimensões lineares 3 vezes maiores? 


P. 2 — Imagine um sistema solar semelhante ao exis- 
tente com dimensões lineares 3 vezes maiores e com 
as massas do Sol e dos planetas 5 vezes menores. Qual 
é a relação entre os periodos de translação dos planetas 
no sistema real e no sistema imaginado? 


P. 3 — Uma ponte tem um vão de 120m. Qual deve 
ser, num planeta em que a aceleração da gravidade seja 
10 vezes superior à da Terra, o vão de uma ponte 
semelhante à anterior, feita com o mesmo material, e tal 
que nela seja a mesma a fadiga máxima devida ao peso 
próprio? 


P. 4 — Considere dois sistemas de condutores no vá- 
cuo, o segundo com dimensões lineares três vezes supe- 
riores às do primeiro. Suponha que as cargas dos con- 
dutores do 2º sistema são sete vezes menores do que 
as cargas dos condutores do 1.º sistema. Determine a 
relação entre as resultantes das forças electrostáticas 
que se exercem sobre os condutores do primeiro e sobre 
os condutores do segundo sistema. Determine, ainda, as 
relações entre os coeficientes de potencial e de capa- 
cidade dos condutores de um e de outro sistema. 


P. 5 — Um iman exerce sobre uma peça em ferro uma 
força de 1 N. Suponha um outro iman e uma outra peça 
em ferro, de formas semelhantes, de dimensões lineares 
duas vezes maiores, feitos de materiais idênticos, e colo- 
cados numa posição semelhante. Que força se exerce 
entre eles? 

Resoluções (abreviadas): 


1 — Além das dimensões geométricas as grandezas 
que interessam à resolução do problema são a densi- 
dade p e o módulo de elasticidade E de dimensões 
(]=M L% e [El=ML! T2 As condições L | =3; 


(12) 
à 2 1 item-l 
M (12) É 1; M (2) Lazy Triay = | permitem-nos ime- 
diatamente obter T,,2, = 3. 


1 
(12) + 
Sendo a dimensão da frequência [f]=T'! temos 
A frequência da vibração da 2.º barra é 24 s'!, 


er 
(12) 


2 — Em astronomia os cálculos baseiam-se na equa- 
ção f=m.a e na lei de Newton. Além das massas e 
das distâncias há a considerar a constante de gravi- 
tação G de dimenões [G] = M:! L? T? que como cons- 
tante universal é a mesma num caso e noutro. Temos 
assim Gy = 1. As condições Lp, =3 May =5"; 
Matas Leça) Tia = 1 dão-nos quase imediatamente 
Taz = (5.3*)'8. 


Pa: a dus = Sd ai 
3—De: 93 =L02» "y= 10. Pao =Maa 3) =! 


== = 2 = = = 101 
e ay Mas Las agp! obtemos 1,,5,=L(,5,=10 
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4— Temos de considerar a lei de Coulomb onde 
figura a constante universal E, de dimensões [E] = 


= M1! LS Tº |*. Para sistematizar os cálculos é útil 
organizar um quadro como o que se segue: 


d L L(42) 

q Et q (123)= 12) Toa) = 7" 

E (MIL TP | eo Mão Elo Tha amy = 
F IML TZ E = Mao Lay Tán=? 

Cj |Mº L2 TP) Cito = Es 
Cii IM L2 TA |-2 | Cia = a 


Às três primeiras linhas permitem-nos obter 


M 15) Lo 2) = 32.7? valor que substituido na 4.b linha nos 
dá f,,2, = 1/32.7º. Dum modo semelhante calculariamos 
Cy =3€ Ciy= 1/3 


5 — Admitimos que a magnetização do iman é dada 
por M=M,+Xn He a da peça em ferro por uma 
fórmula semelhante mas sem magnetização permanente, 
As forças dependem da geometria do sistema e das 
grandezas M,, ti, e X- Para resolver o problema pode- 
mos fazer um quadro semelhante ao anterior. 


t L (12) = 

ML T+? |-42 Boy= M, 12) La) Tê rá = 
1 
Lt) Moy= (42) 42) = 

M LT? fo) = M, 12) Eq) Ta = 2º 


As três primeiras linhas permitem-nos obter 
Mais Lad Lis = 2? As forças quadruplicam quando 
as dimensões lineares duplicam. 


Detectadas as grandezas que interessam à resolu- 
ção do problema, obtemos o resultado quase automati- 
camente. Notemos que encontramos neste problema 
Ho) = 1. No caso da magnetização ser devida também 
a correntes e dos materiais obedecerem a leis M=M (H) 
não lineares já nada poderiamos fazer, Deixava de haver 
semelhança entre os dois sistemas, embora feitos de 
peças geometricamente semelhantes e dos mesmos ma- 
teriais. 


8 — Um problema particular 


O problema que a seguir apresentamos leva-nos a 
considerar uma noção de semelhança, ligeiramente mo- 
dificada, utilizável nalguns casos especiais. 


— Um ponto de uma estrutura desloca-se de 7 mm 
quando ela é carregada com determinado conjunto de 
cargas. 
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Imagine uma estrutura semelhante à anterior, feita 
com os mesmos materiais, e de dimensões lineares 
duas vezes maiores. 


Quanto se desloca o ponto da segunda estrutura, 
homólogo do ponto da primeira, quando ela é carregada 
com um conjunto de cargas, semelhantes às da primeira 
estrutura, mas de intensidades três vezes maiores? 


Resolução: 


Representando por d as dimensões lineares das 
estruturas e por 8 os deslocamentos dos seus pontos 
devidos às cargas, temos, num quadro semelhante aos 
dos problemas anteriores; 


dz) = (12) ad é 


õ L dy= Lay=? 

E | 

[o] ML! T2 Sao Mp) La) Tão = Ê 
E M L! T Emos M (12) La) Ty= 
f ME “Te Fa» May & Tág=3 


mm) (2) 


Para haver semelhança (no sentido da definição 
atrás adoptada) teriamos de ter na segunda linha: 


d(»)= Ly=2 
o que implicaria na última linha: 


Micas bzgãs Tm 


Fa = 2 “qa qa 


o que não é concordante com um dos dados do pro- 
blema. 


Os dois sistemas físicos considerados (as duas 
estruturas com as respectivas cargas) não são pois sis- 
temas semelhantes. 

Mas o problema não deixa de ser resolúvel, 

Sabemos que numa estrutura (na hipótese linear) 
os deslocamentos são proporcionais às cargas. Os sis- 
temas físicos SF, formados por uma dada estrutura com 
diferentes hipóteses | de carga semelhantes entre si, 
são sistemas físicos não semelhantes dado que por um 
lado temos sempre d,,;, = L,yy = 1 e por outro 8,4 = 
= Ls = fá, 51. Não deixa, no entanto, de haver 
semelhança entre os deslocamentos dos diferentes sis- 
temas e entre as forças dos diferentes sistemas. 

A condição quebrada é a da igualdade entre os 
coeficientes de proporcionalidade (na passagem de um 
sistema para o outro) correspondentes às distâncias glo- 
bais e aos deslocamentos, grandezas com a mesma 
dimensão fisica. 

Tudo se passa como se houvesse semelhança, mas 
estas duas grandezas tivessem dimensões diferentes. 

Escrevendo [d]= L e [5] = L*, as fórmulas da Me- 
cânica e da Elasticidade permitem-nos escrever as di- 
mensões das diferentes grandezas que nos interessam 
num sistema de grandezas fundamentais M, L, T e L*. 
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Usando estas fórmulas, e seguindo a técnica dos 
problemas anteriores escrevemos: 


d L day=Lç2=2 

ô L.* dy =Lfh2 e 

E a Si Ea) Lay Lo =? 

oi Mt 4 "(o = Ma Lia) Ta» cilia 
E MiLta Ts ES M (12) La) Tãa = l 
f M LT? Fay M (42) L 42) Taz =o 


Os dados da 1.º, 5.º 6.º linhas permitem-nos obter 
Lfs2) L (12) = 3, (= L* == 3ja, Eq =% . Temos, final- 
mente 8 (12) = 105 mm. 


Na resolução deste problema fizemos uso de uma 
noção de semelhança parcial. Entre dois sistemas há 
semelhança parcial quando (satisfeitas todas as outras 
condições de semelhança) duas grandezas de igual di- 
mensão não satisfazem à condição 3 da secção 2, por 
terem coficientes de proporcionalidade diferentes na 
passagem de um sistema para outro, sendo, no entanto 
satisfeitas as condições 3 quando as duas grandezas 
são consideradas de dimensões independentes. 


Il — O TEOREMA «Pl» OU DE BUCKINGHAM 
1 — Comentário prévio 


Consideremos a seguinte afirmação que designare- 
mos por afirmação À: 


«As equações da Fisica são invariantes numa mu- 
dança de sistema de unidades.» 


Que significa esta afirmação? 

Conhecidas as fórmulas de transformação dos va- 
lores numéricos dum conjunto de grandezas g!,...,. g* 
numa mudança de sistema de unidades, a passagem 
da uma relação do tipo 


1 k 
Hlg; ...g,)=0 


entre os valores numéricos das grandezas num sistema 
de unidades 1 para a relação correspondente 


k 
lg... 92) =0 


entre os valores numéricos das grandezas num sistema 
de unidades 2 é um elementar problema de matemática. 

A afirmação A diz pura e simplesmente que expri- 
mindo-se uma lei da Fisica por equações como as acima 
indicadas, as funções f, e f, são iguais. 

Vimos que os sistemas de unidades eram exacta- 
mente escolhidos de modo a que tal fosse verdade no 
caso das leis que designamos por leis iniciais da Fisica. 
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No caso da gravitação vimos que a expressão da 
força em função das massas e das distâncias não era 
invariante mas que a invariância era obtida por meio 
d.: introdução de uma constante universal, 

Em geral, uma lei não invariante que num dado sis- 
tema de unidades se exprime por uma equação do tipo 


1 k 
fila: assa 91) = 0 
pode ser sempre transformada numa lei invariante por 
meio da introdução de constantes c! que não necessitam 
de ser em número superior a k porque sendo c' de 
dimensão igual à grandeza g' e de valor cj= 1 no sis- 


tema adoptado. 
1 | 
' 91 gr 
9) = joe —— E= O) 
Cc, Cc: 


já é uma lei invariante, (No caso de algumas grande- 
zas g! terem a mesma dimensão as respectivas cons- 
tantes ci podem identificar-se). 

Apresentadas assim as coisas a afirmação A apare- 
ce-nos como inteiramente vazia de conteúdo físico. 

Na realidade, no estudo dos fenómenos fundamen- 
tais o número de constantes c' necessárias para obter a 
invariância das leis é sempre inferior aos das grande- 
zas de dimensão diferente em jogo. A invariância obtém- 
-se pois sem necessidade de fazer uso do número má- 
ximo de constantes em principio eventualmente necessá- 
rias, 

É este facto subjacente que confere um conteúdo 
físico (importantissimo) à afirmação A. 

É fácil ver, tendo em conta o que atrás dissemos 
sobre semelhança, que no caso extremo de ser necessá- 
rio introduzir o número máximo de constante ct para 
obter a invariância de uma lei, não é possivel, no domi- 
nio de aplicação dessa lei, reproduzir situações seme- 
lhantes a escalas diferentes. 

É o que se passa com frequência no estudo de 
alguns fenómenos particulares, caso em que os ck não 
são constantes universais, mas sim constantes caracte- 
risticas dos materiais. Basta, por exemplo, no estudo 
de um fenômeno magnético ter um material em que a 
relação B = f(H) é não linear, para não ser possivel 
reproduzir campos magnéticos semelhantes mas diferen- 
tes (1H (12) E 1). 

E, no próprio caso dos fenómenos fundamentais, 
podemos suspeitar que a semelhança só se verifica nos 
limites de uma escala muito larga no interior da qual 
o efeito de parte das constantes c* não se faz sentir. 
Nunca haveria assim, em rigor absoluto, situações seme- 
lhantes a escalas diferentes. 


f(c'!, ' cr - 9, 1 


2 — Demonstração 


Partamos agora da afirmação A cujo significado foi 
esclarecido. 

Admitamos que sabemos que numa certa situação 
fisica o valor de uma grandeza g! é determinado pelos 
valores de n-1 grandezas g?,....g". Representemos por 
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Gi os valores destas grandezas num dado sistema de 
unidades. O nosso objectivo é conhecer a relação 


Lj 


1 E. 
G, = h (G, Esp G,) 


a que podemos dar a forma: 


Sendo as leis da Fisica invariantes (o que significa 
que nelas incluímos todas as constantes necessárias 
para assegurar essa invariância) entre os valores G; 
das grandezas q' num segundo sistema de unidades 
deve verificar-se a mesma relação 


RG; ..., 65) =0 


Entre as n grandezas g”" podemos sempre escolher 
um conjunto de grandezas que podem ser tomadas como 
grandezas fundamentais, e em relação às quais as res- 
tantes se podem exprimir dimensionalmente, Admitamos 
que essas grandezas são, por exemplo, as três primei- 
ras. Representemos por 4, B e C as suas dimensões. 
Podemos escrever para todas as outras grandezas as 
equações dimensionais 


k ak pk 
lgJ=A B CC 


em que os expoentes são faceis de calcular. Conside- 
remos os números 


k G 
uh | ak é Pk vk 
(G,) (GG) (Gi) 
(os três primeiros manifestamente iguais à unidade) 
independentes do sistema de unidades inicialmente esco- 
lhido. 

É fácil ver que, num sistema de unidades em que 
os valores correspondentes às três primeiras grandezas 
do sistema em estudo sejam a unidade as restantes 
grandezas tomam os valores II*. Como a invariância 
se verifica ainda, necessariamente, para este novo sis- 
tema é válida a relação 


1, 14 0, 0)= FOR...) = 0 


Vemos que o problema da procura de uma relação 
funcional entre n variáveis se reduz ao problema da 
procura de uma relação entre n-s variáveis (s — número 
de grandezas que podem ser tomadas como fundamen- 
tais entre as grandezas em jogo). Este resultado é o cha- 
mado teorema «Pi» ou teorema de Buckingham, 

Nas exposições correntes, este teorema é demons- 
trado a partir da afirmação A, em geral aceite à partida 
sem qualquer discussão prévia, Na falta desta discussão 
o teorema aparece como um mero resultado de mate- 
mática cujo significado físico nos escapa. 

Em particular é de sublinhar que os números II 
correspondentes a situações semelhantes (quando estas 
existem) são iguais. Mas este é, de certo modo, um 
resultado final e restrito, e não nos parece que a técnica 
dos números «Pi» seja a técnica mais conveniente para 
abordar em toda a sua globalidade o problema da seme- 
lhança. 


(Continuação da pág. 26) 
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LUÍS VALADARES TAVARES — Lima perspectiva 
sobre os problemas pedagógicos no IST, «Técnicas, 
Lisboa, LV (460), Outubro 1980, p. 5-8. 


JOSÉ CARLOS DE O. S. HORTA — As crostas 
calcárias e a construção de pavimentos de estradas 
na África do Norte [primeira parte), «Técnica», Lis- 
boa, LV (460), Outubro 1980, p. 9-22. 


Em clima mediterrânico sub-húmido a semi-árido 
a calcite é dissolvida nos maciços montanhosos pelas 
precipitações da estação fria e deposita-se, no verão, 
nos solos das planícies, cujos perfis se diferenciam à 
medida que se vão enriquecendo em calcário, À, RUEL- 
LAN deu, em 1970, uma descrição genética e uma 
classificação muito completas dos diferentes estados de 
acumulação de calcário nos solos, As crostas calcárias 
folheadas, podendo compreender uma lage petrificada 
na parte superior, correspondem ao que se tem cha- 
mado tradicionalmente tufo na África do Norte. Esta 
designação incorrecta foi, aqui, consagrada pelo uso em 
técnica de estradas, Trata-se, com efeito, dum exce- 
lente material de pavimentação quando utilizado den- 
tro das condições climáticas que correspondem ao en- 
crostamento dos perfis pedológicos. 


FERNANDO ALBERTO GONÇALVES — Sobre a re- 
presentação analítica de resultados experimentais, 
«Técnica», Lisboa, LV (460), Outubro 1980, p. 
23-26. 


À linearização de relações funcionais permite a fácil 
representação analítica de resultados experimentais me- 
diante o cálculo de apenas dois parâmetros. 

Quando a relação funcional não é linearizável, por 
exemplo, por se desconhecer a sua forma, e os resulta- 
dos experimentais em coordenadas rectangulares suge- 
rem uma ligeira curvatura, a técnica de abaixamento 
de grau permite dar forma linear a uma relação qua- 
drática, o que simplifica grandemente os cálculos dos 
coeficientes da equação procurada. 


Ano LV — Outubro 1980 


J. J. DELGADO DOMINGOS, M. FÁTIMA PINTO 
e M. TERESA PONTES — Ocorrência média anual 
no território português das classes de estabilidade 
atmosférica Pasquill-Giftord, «Técnica», Lisboa, LV 
(460), Outubro 1980, p. 27-42, 


Com o objectivo de fundamentar a regulamentação 
aplicável às actividades que contribuem para a poluição 
atmosférica no Território Nacional, efectuou-se a carac- 
terização da ocorrência média anual das classes de 
estabilidade atmosférica Pasquill-Giflord no Território 
Continental. 

Esta caracterização permitirá estabelecer as expres- 
sões mais adequadas para o cálculo da dispersão de 
efluentes atmosféricos no Pais. 


A. BETÂMIO DE ALMEIDA — Regime variável em 
pressão: Notícia sobre bibliografia recente, «Téc- 
nica», Lisboa, LV (460), Outubro 1980, p. 43-44, 


ANTÓNIO BROTAS ——- Notas sobre a teoria da se- 
melhança, «Técnica», Lisboa, LV (460), Outubro 
1980, p. 45. 


Após discussão do problema da invariância das 
equações da Física nas mudanças de sistemas de uni- 
pasro é apresentada uma definição de semelhança 
SICa. 

É em seguida apresentado o problema das possibi- 
lidades de realização de sistemas fisicos semelhantes, 

Uma notação proposta neste texto permite a reso- 
lução quase automática dos problemas de semelhança 
(das relações entre as - prandezas físicas de sistemas 
semelhantes). 

Estas notas foram redigidas para serem usadas pe- 
los alunos da cadeira de Introducão à Física dos Meios 
Contínuos do 1, 8. T. Com a sua publicação na «Téc- 
nica» pretende-se contribuir para que nas cadeiras de 
Fisica das escolas de Engenharia passe a zer dada 
peça atenção a este assunto e a este tipo de pro- 

EMAS. 
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Synopsis of articles published in «Técnica» n.º 460 
LV — October 1980 


J. J. DELGADO DOMINGOS, M. FÁTIMA PINTO 
e M. TERESA PONTES — Annual mean occurence 
of Pasquill-Gifford atmospheric stability classes in 
Portugal, «Técnica», Lisboa, LV (460), Outubro 
1980, p. 27-42, 


Annual mean oceurence of Pasquill-Gilford stability 
classes in Portugal was made in order to fundament 
the policy to rule atmospheric pollutant activities. 

Based on this caracterization more correct expres- 
sions to calculate the dispersion of atmospherie pullu- 
tants in Portugal can be set. 


A. BETÂMIO DE ALMEIDA — Unsteady pressure 
flow: Recent books review, «Técnica», Lisboa, LV 
(460), Outubro 1980, p. 43-44. 


ANTÓNIO BROTAS — Notes on similitude theory, 
«Técnicas», Lisboa, LV (460), Outubro 1980, p. 45. 


- After discution of the problem of the invariance of 
Physical equations in the changes systems of unit, a 
difinition of physical similarity is presented. 

Then the problem of possibility of realization of 
similar physical systems is presented, 

- À notation proposed in this text gives the possi- 
bility of an almost automatic resolution of the pro- 
blems of similarity (of the relations between physical 
variables of similar systems). 

These notes were written to be used by scholars 
of «Introduction to the Physics of Continous Medias 
at the TI. 8. T. With the notes here presented we are 
trying to given a contribution so that in several en- 
gineering schools a greater attention to this subject 
and this kind of problems shall be given. 
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LUÍS VALADARES TAVARES — 4 discussion on 
pedagogic issues for [S.T., «Técnica», Lisboa, LV 
(460), Outubro 1980, p. 5-8. 


JOSÉ CARLOS DE O. S, HORTA — Calcareous 
crusts (caliches or calcretes) and road pavement 
construction in North Africa fpart one), «Técnica», 
Lisboa, LV (460), Outubro 1980, p. 9-22. 


Under the subhumid and semi-arid mediterranean 
climates calcite is solubilized on the mountainous 
massifs by the cold season precipitations to be depo- 
sited in summer in the soils of the plains, which profi- 
les differentiate as their enrichment in limestone 
proceeds, A. RUELLAN has given in 1970 a complete 
genetic description and classification of the different 
levels of this accumulation of limestone in the soils. 
Foliated caleretes with eventually a petrefied hardpan 
at their upper part correspond to the materials which 
have been traditionnally designated by tufa in North 
Africa. This designation is not correct but has been 
imposed by routine in road technology. It is, indeed, 
an excellent pavement material whether used within 
the climate conditions corresponding to the encrusting 
of the pedological profiles, 


FERNANDO ALBERTO GONÇALVES — On the ana- 
lytical description of experimental results, «Técni- 
ca», Lisboa, LV (460), Outubro 1980, p. 23-26. 


The linearization of functional relations allow an 
easy representation of experimental results through the 
calculation of only two parameters. 

When the functional relation is not linearizable, for 
example if we do not know its form, and the results 
in rectangular coordinats suggest a slight curvature, a 
lowering of degree permit to linearize a quadratic rela- 
tion, which simplifies greatly the calculation of the 
parameters of the searched equation, 
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PEDIDOS A [BENS 


If you could talk to Orville Wright, he'd 
tell you the problems he faced as a turn- 
of-the-century engineer. You could tell him 
all about the technological solutions 


available to today's engineer and scientist... 


particularly a 20th century phenomenon 
that tests assumptions and defines models 
before a project gets off the ground. 

The Apple personal computer. 


Computation, calculation, 
analysis...the power to 
pilot your projects. 


With a highly-integrated system from 
the extensive Apple personal computer 
family, Orville and brother Wilbur would 
have increased their productivity. Perhaps 
even launched the Kitty Hawk Flyer well 
before 1903. 

An Apple in their hangar would have 
freed them from the time and tedium of 
crunching numbers by hand. 

An Apple in your lab or office will give 
you the problem-solving capabilities you 
demand from a big computer...without 
the time-consunung problems typical of 
remote processing. 

But the Apple system solution doesn't 
stop there. Itkeeps on soaring with proven 
performance, power and expandability 


Apples existing software library includes 
a program that plots the shape of an 
airfoil, given its parameters. 


that's unparalleled for analyzing alterna- 
tive paths of design and modeling a 
wide variety of physical processes. 

Want more memory? Depending on 
your choice of system, Apple has memory 
expandable to 64K bytes or 128K bytes. 
Prefer wide displays? Choose 40 or 80 
characters. Need to control instruments 
in the lab? Get on the IEEE 488 bus. Over 


100 companies also supply peripherals for 
Apple because Apple is the most popular 
personal computer with the least com- 
plicated interface. 

Want an efficient system of data storage 
and access? Apple's 5 /4" disk drive not 
only offers you increased application 
versatility, but high density (143K bytes), 
high speed and low cost. You can even add 
up to four or more drives to your 
Apple system. With 
proven reliability; 
no wonder it's the 
most popular drive 
on the market today. 
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Wilbur determined that birds didn't have to constanti y flap 
their wings to fly. With an Apple, be could've determined the 
fixed-wing design of the Kitty Hawk Flyer much faster 


FORTRAN that helped 
to design a 20th century 
flying machine. 

Fluent in the same language that helped 
to design the 747, Apple FORTRAN lets you 
tackle differential equations at the touch of 
a key. And since more the 170 companies 
also offer software for the Apple family, 
you can have one of the most impressive 
program libraries ever...including vast 
subroutine libraries for math, science, 
engineering and statistics. When you write 


o 
“PP Ãovourown 
” programs, the 
Apple also speaks 
in languages other 
than FORTRAN: Pascal, 
BASIC, PILOT and 6502 

assembly language. 


Where to learn 


js 


the small-yet- 
serious 
solution. 
Let your 
imagination 
soar wilh 
Apple. 
| Discover the 20th 
century tool versatile 
enough to monitor 
quality controls and 
manufacturing sched- 
ules, orchestrate 
tolerance tests and 
determine alternative 
parts selection. Learn why Apple emerges 
as the technological leader of reliable 


personal computer products that increase 


your productivity 


Visit our sole distributor 

in Portugal, Sorubal S.A.R.L,, 
Rua Gen. Pimenta de Castro 
15 8.º, 1700 Lisboa, 

or call 896555 (Lisbon) 
telex 12775 SORBAL P, 


more about Apple, 


